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AERODYNAMIC  DRAG  OF  CYLINDRICAL  VEHICLES  MOVING 
CONCENTRICALLY  WITHIN  LONG  TUBES 


1.  INTRODUCTION 


The  objective  of  this  study  is  to  model  the  flow  field  in  the  annular 
region  between  a  moving  cylindrical  vehicle  and  a  concentric  tube.  The 
results  of  the  analysis  may  be  applied  to  both  the  swimout  and  impulse 
launch  methods  (figure  1-1).  These  methods  are  typically  used  by  the  Navy 
in  launching  various  projectiles.  The  figure  shows  the  associated  velocity 
profiles  in  the  annulus  for  each  launch  method.  For  swimout  launch,  the 
direction  of  fluid  motion  is  opposite  that  of  the  vehicle,  which  gives  rise 
to  shear  stresses  that  oppose  the  motion  of  the  vehicle.  During  impulse 
launch,  however,  the  annular  profile  shows  that  the  fluid  moves  in  the  same 
direction  as  the  vehicle,  and  the  shear  stresses  aid  launch. 

Figure  1-2  shows  the  cylindrical  vehicle  during  swimout,  the  vehicle 
having  moved  from  the  initial  position  x  to  some  new  position,  x  +  dx.  The 
vehicle  displaces  a  quantity  of  fluid  that  must  accelerate  through  the 
annular  gap  to  fill  the  void  created  at  the  rear  of  the  vehicle.  The  shear 
stresses  resist  vehicle  motion  and  become  more  pronounced  with  increasing 
vehicle  speeds  and  decreasing  annular  gaps. 

For  the  alternate  scenario  of  impulse  launch,  the  fluid  is  forced 
through  the  annulus  from  the  rear  of  the  cylinder.  The  shear  stresses  pull 
the  vehicle  along  the  direction  of  flow.  Since  a  high-pressure  region  exists 
behind  the  cylinder  relative  to  the  forward  end,  a  force  imbalance  exists 
that  consequently  induces  vehicle  movement.  As  gap  clearances  become  small, 
the  following  events  occur  that  aid  the  impulse  launch: 

1.  Annular  fluid  velocities  increase,  magnifying  the  wall  shear 


stresses. 

2.  The  volume  flux  of  blowby  decreases,  resulting  in  a  larger  pressure 
differential  across  the  vehicle. 


1-1 


Experimental  efforts  to  understand  turbulent  annular  tube  flow  were 

reported  as  early  as  1907  (reference  1).  In  1968,  reference  2  documented 

the  successful  prediction  of  friction  factors  for  fully  developed,  turbulent 

annular  tube  flow.  The  theoretical  results  presented  were  for  the  range 

4 

6000  <  Re^  <  4.5  x  10  ,  for  six  radius  ratios  ranging  from  1.05  to  50.* 
Deissler's  equation  was  used  to  model  the  eddy  diffusivity  and  obtain 
velocity  profiles  for  the  sublayer,  while  Von  {Carman's  similarity  hypothesis 
was  used  to  obtain  profiles  in  the  turbulent  layer.  The  combined  use  of 
these  two  equations  resulted  in  reasonable  agreement  with  the  experimental 
results  previously  obtained  by  Quarmby  in  reference  3.  Then,  in  1971,  a 
model  was  developed  for  the  hydrodynamic  boundary  layer  growth  in  the  entry 
region  for  annular  tube  flow  using  the  momentum  integral  technique 
(reference  4).  The  velocity  profiles  were  determined  by  the  use  of 
Reichardt's  expression  for  eddy  diffusivity,  and  were  adjusted  for 
calculations  near  the  wall  by  implementation  of  Van  Drieat's  damping  factor. 

In  1981,  the  successful  analytical  models  of  references  2  and  4  were 
combined  and  extended  to  include  a  moving  boundary  (that  is,  a  venicle 
moving  within  a  concentric  tube)  as  reported  in  reference  5.  The  predicted 
vehicle  drag  coefficient  was  16  percent  lower  than  the  experimental  results 
of  previous  investigations,  as  compiled  by  Davidson  in  reference  6.  Davidson 
reduced  and  correlated  these  results  to  a  common  plot  of  drag  coefficient 
based  on  vehicle  speed  versus  the  area  blockage  ratio.  In  reference  5,  Sud 
and  Chaddock  stated  that  neglect  of  entrance  effects  and  the  use  of  a 
simplified  model  for  the  wake  region  in  their  analysis  may  have  contributed 
to  the  discrepancy  between  their  results  and  the  correlation  of  Davidson. 

The  calculations  of  Sud  and  Chaddock  were  limited  to  a  single  tube/ 
vehicle  radius  ratio  for  two  Reynolds  numbers.  The  analytical  approach  and 
general  iterative  scheme  of  this  present  study  are  nearly  identical  to  those 
used  by  Sud  and  Chaddock.  In  the  present  effort,  however,  two  major  improve¬ 
ments  were  implemented  that  have  been  developed  through  exhaustive  numerical 
experimentation.  These  improvements  consist  of  curve  fits  to  estimate  the 

*Symbols  are  defined  in  the  "List  of  Symbols,"  p.  vi. 


initial  parameter  values  and  a  series  of  weighted  linear  interpolations  to 

steer  the  iteration  process.  These  techniques  enable  convergence  to  a 

4  7 

solution  for  a  wide  range  of  Reynolds  numbers*  10  <  Rey  <  10  ,  and  radius 

ratios*  1.01  <  b  <  2.0*  while  keeping  the  volume  flux  constant  within  0.01 
percent.  Additional  techniques  were  developed  for  solving  the  developing 
flow  problem.  These  are  discussed  under  "Analytical  Survey"  in  section  3  of 
this  report. 
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2.  ANALYSIS  OF  FULLY  DEVELOPED  FLOW 


EXPERIMENTAL  SURVEY 


Steady,  fully  developed  annular  tube  flow  was  experimentally 
investigated  as  early  as  1907  (reference  1).  Since  that  time,  many 
investigators  have  conducted  experiments  to  characterize  this  flow  because 
it  occurs  frequently  in  heat  exchangers.  Of  those  investigators,  some  found 
that  r^/r^  >  l,  whereas  others  found  that  <  1*  This  particular 

characteristic  is  perhaps  the  most  controversial.  These  quantities  are 
depicted  in  figure  2-1. 

More  recent  results  (reference  7)  confirm  those  investigators  who 
found  K  l»  which  now  seems  to  be  the  generally  accepted  result. 

Consequently,  only  references  3,  8,  9,  and  10  will  be  discussed  in  this 
survey . 

In  reference  8,  the  flow  of  air  was  measured  through  two  annuli,  namely 

b  *  6.17  and  1.54  for  4000  <  Re^  <  22,000.  The  authors  found  r^/r^  »  1 

for  all  but  Re^  *  1820  with  b  *  6.17.  Velocities  at  points  beyond  rmt 

(near  the  tube)  were  found  to  agree  with  the  laminar  correlations  of 

reference  11,  unlike  velocities  below  r  ... 

mt 

In  reference  9,  flow  of  air  was  measured  through  annuli  in  which  b 
ranged  from  2  to  5.2  for  10,000  <  Re^  <  150,000.  The  investigators  arrived 
at  the  following  correlation  for  r  . : 

DC 


b.  TURBULCNT  FLOW 


Figure  2-1.  Deplot  Ion  of  rat,  raL  for  Fully  Developed  Flow  Through 

Conoentrlo  Annuli 


The  experimenters  in  reference  10  used  air  to  test  the  flow  through 
annular  pipes,  with  7  radius  ratios  in  the  range  of  1.78  <  b  <  16  for  the 
range  47,000  <  Re^  <  327,000.  Friction  factors  were  found  to  be 
between  1  and  10  percent  above  pipe  flow  data,  with  little  dependence  on  b. 
The  authors  concluded  that  friction  factors  are  independent  of  Reynolds 
numbers  for  the  range  of  radius  ratios  and  Reynolds  numbers  tested. 

According  to  reference  9,  the  experiments  reported  in  references  12  and 
13  were  performed  with  great  care.  Reference  12  found  rm^rnj  <  1 
Reynolds  numbers  around  250,000  with  b  as  large  as  19.  The  working  fluid 
was  water.  Reference  13  presented  data  for  Reynolds  numbers  from  10,000  to 
70,000  for  air,  and  for  the  radius  ratio  of  1.845. 

Reference  3  tested  three  radius  ratios:  2.88,  5.62,  and  9.37  for  a 

Reynolds  number  range  from  6000  to  90,000,  using  air  as  the  working  fluid. 

The  author  found  the  difference  between  r  .  and  r  .  to  be  a  function  of  both 

mt  ml 

Reynolds  number  and  radius  ratios  for  40,000  <  Re.  <  50,000,  and  to  be  a 

o 

function  of  b  only  for  larger  values  of  Re^.  Friction  factors  were  found 
to  be  nearly  independent  of  radius  ratios,  and  were  closer  to  plain  tube 
experimental  results  than  those  of  reference  10. 

In  reference  7,  three  annuli  with  radius  ratios  of  2.53,  5.68,  and 
11.36  were  tested  for  the  Reynolds  number  range  from  20,000  to  300,000. 
Friction  factors  were  found  to  be  8.5  percent  above  pipe  flow  for  b  *  11.36 
and  only  5  percent  higher  for  the  ratios  3.53  and  5.68.  It  was  also  found 
that  a  considerable  portion  of  the  boundary  layer  associated  with  the  outer 
wall  could  be  modeled  using  the  law-of-the-wall.  The  profiles  within  the 
inner  boundary  layer  could  be  described  using  modified  law-of-the-wall 
relations  for  b  ■  2.53  and  5.68.  However,  for  b  »  11.36,  law-of-the-wall 
relations  could  not  be  deduced. 


ANALYTICAL  SURVEY 

In  reference  2,  Quarmby  disagreed  with  Tomotika  and  Imai  (reference  14), 
who  found  >  And  stated  further  that  all  other  investigators  known 
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to  him  found  r^/r^  <  1,  with  the  exception  of  those  in  reference  8,  who 
found  the  ratio  rmt/ rm^  t0  b«  either  less  than  or  equal  to  unity  (for  low 
Re^).  Quarmby  used  the  Von  Kerman  similarity  hypothesis  and  Deissler's 
expression  for  eddy  diffusivity  to  obtain  turbulent  layer  and  sublayer 
velocity  profiles,  respectively.  His  analytical  results  compared  well 
with  his  experimental  results  in  reference  3  for  6000  <  Re^  <  45,000,  with 
1.05  <  b  <  50. 

In  reference  15,  Quarmby  further  investigated  the  Re^  effect  on 
u+(y+)  profiles  for  the  limiting  cases  of  parallel  plate  and  circular 
tube  flows.  He  corrected  deficiencies  in  his  earlier  analytical  model 
(reference  2)  where  the  shear  stress  was  assumed  uniform  across  the  duct, 
which  is  analogous  to  neglecting  the  molecular  viscosity.  Further,  the 
condition  of  a  zero  profile  gradient  at  r  previously  unfulfilled  in 
reference  2  was  also  corrected.  When  modifications  of  reference  15  were 
applied  to  Deissler's  and  Von  Karman's  equations  for  eddy  diffusivity, 
good  agreement  was  found  with  the  experimental  results  of  reference  16  for 
the  limiting  cases  of  plain  tube  and  parallel  plate  channel  flows. 

Up  to  this  point  of  the  survey,  no  provisions  had  been  made  by  any 
investigator  to  include  the  effects  of  a  moving  wall  in  annular  tube 
flows.  Sud  and  Chaddock  (reference  5)  used  Quarmby' s  modification  in  the 
sublayer  (reference  15)  and  the  turbulent  layer  expressions  of  reference  2 
in  the  extended  problem  of  a  vehicle  moving  concentrically  within  a  tube. 
Good  agreement  was  found  with  reference  3  for  the  generic  case  of  no 
moving  boundaries.  An  mentioned  in  section  1,  reference  5  also  compared 
well  with  reference  6.  In  view  of  the  success  of  Quarmby  and  Sud  and 
Chaddock,  and  the  fact  that  their  analyses  were  physically  the  most 
fundamental,  an  approach  similar  to  that  of  reference  5  was  used  in  this 
study. 

ASSUMPTIONS  AND  BOUNDARY  CONDITIONS 

In  this  study,  the  vehicle  was  assumed  to  move  within  the  concentric 
tube  at  a  constant  velocity.  As  in  reference  5,  the  analysis  was  simplified 


by  attaching  the  coordinate  system  to  the  vehicle  surface  just  aft  of  the 
nose  section,  thus  making  the  flow  field  within  the  annulus  appear  steady 
(figure  2-2).  If  the  coordinate  system  were  attached  to  the  tube,  the  flow 
field  observed  would  change  with  the  vehicle's  passage,  and  the  following 
sequence  of  events  would  occur.  As  the  vehicle  approaches  the  coordinate 
location,  a  pressure  rise  would  occur.  During  the  vehicle's  passage, 
annular  effects  would  be  present.  Finally,  as  the  rear  of  the  vehicle 
crosses  the  coordinate  location,  the  pressure  would  drop  rapidly,  leading 
into  wake  effects.  These  events  would  be  time' dependent  and  would  add 
unnecessary  complications  to  the  analysis  if  included  as  a  result  of 
coordinate  placement.  Because  the  coordinate  system  was  attached  to  the 
vehicle  surface,  the  tube  appears  to  be  moving.  Figure  2-3  shows  how  the 
velocity  profile  would  be  interpreted  by  virtue  of  coordinate  placement  on 
either  tube  or  vehicle  wall.  Note  that  the  shapes  of  these  velocity 
profiles  are  in  fact  identical,  with  the  exception  of  the  position  of  the 
zero  velocity  vector  and  consequently  the  total  volume  flux. 

The  velocity  profile  at  the  annular  entrance  was  assumed  uniform.  This 
is  a  reasonable  assumption  provided  the  boundary  layers  are  initially  tur¬ 
bulent  and  the  vehicle  nose  is  well  rounded,  as  discussed  in  reference  4. 
Entrance  profiles  are  actually  nonuniform  because  of  boundary  layer  growth 
on  the  vehicle  nose  section  and  may  result  in  larger  drag  estimates  if 
properly  accounted  for.  This  is  especially  true  for  blunt  nose  vehicles 
where  separation  is  likely  to  occur.  The  following  assumptions  were  made: 

1.  Steady  flow  is  present,  i.e.,  dV/dt  »  0. 

2.  There  is  a  uniform  velocity  profile  at  the  annular  entrance.  This 
is  a  reasonable  assumption  provided  a  rounded  nose  section  is  being  modeled 
and  the  boundary  layers  are  turbulent  from  the  start  (reference  4).  At 

x  •  0.0,  u  ■  Uin  for  all  y. 

3.  Annular  fluid  is  incompressible  with  constant  p. 

4.  There  is  no  blowing  or  suction  at  the  walls.  Therefore, 
u(r)  only...(  no  blowing  or  suction  v  ■  0.0  ). 

3.  Radial  pressure  variations  are  negligible. 

6.  The  tube  moves  with  constant  velocity. 
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The  boundary  conditions  applied  are: 

du/dy)Q  *  0.0. 

0.0  (stationary  vehicle), 
(moving  tube). 


Boundary  condition  no.  1  is  consistent  with  physical  intuition.  The  slopes 
of  the  velocity  profile  for  region  "i"  and  region  "o"  must  be  the  same  at 
r  a  r  with  no  sharp  gradients  or  discontinuities;  furthermore,  the 
slopes  of  the  velocity  profiles  at  this  location  must  be  zero. 


DETERMINATION  OF  VELOCITY  PROFILES 

Figure  2-2  shows  an  enlarged  view  of  the  annular  region.  The  portion 
of  the  flow  field  from  where  the  boundary  layers  start  to  the  point  where 
they  meet  is  termed  the  entry  region  or  developing  region.  Velocity  profiles 
are  functions  of  axial  location  from  the  throat  up  to  this  meeting  point. 
Downstream  of  the  merged  boundary  layers,  the  velocity  profiles  no  longer 
change  with  axial  location.  This  region  is  termed  fully  developed.  The 
fluid  that  enters  the  annular  core  accelerates  up  to  the  end  of  the  entry 
region.  Once  in  the  fully  developed  region,  annular  velocities  have  reached 
their  largest  values,  with  a  maximum  value  of  any  profile  occurring  at 


r  *  r 


mt 


The  analysis  was  dependent  on  determining  the  velocity  profiles,  which 
were  generated  by  implementing  a  two-layer  model.  Deissler's  equation  for 
eddy  diffusivity  of  momentum  was  used  to  obtain  profiles  in  the  sublayer 
where  viscous  shear  is  dominant.  The  eddy  diffusivity  in  the  turbulent 
layer,  which  extends  beyond  the  sublayer  to  the  point  of  maximum  velocity  at 
r  •  r  ,  was  modeled  by  employing  Von  {Carman's  similarity  hypothesis.  As 
previously  mentioned,  this  two-layer  approach  has  produced  friction  factors 
that  are  in  good  agreement  with  experiments. 


2-7 


In  turbulent  theory,  properties  are  generally  represented  in  terms  of  a 
time-averaged  component  and  a  fluctuating  component,  such  as  u  »  u  +  u' , 
where  the  time-averaged  quantity  u  is  defined  by 

T 

u  dt  . 
o 

The  quantity  T  is  the  integration  interval  usually  chosen  to  be  larger  than 
any  significant  time  period  of  the  fluctuations  u*.  In  this  study,  the 
fluctuating  components  were  not  included  so  that  a  mathematical  solution 
would  be  possible.  Average  values  generally  yield  sufficiently  accurate 
results  (reference  17).  Throughout  the  analysis,  all  properties  will 
actually  be  their  averaged  components. 

In  a  one-dimensional  turbulent  flow  field,  the  total  apparent  shear 
stress  (molecular  +  turbulent)  may  be  represented  by 


*  *  0  <v  ♦  '.>  §  •  <2-l) 

where  v  »  e  in  the  sublayer  close  to  the  wall,  and  e  »  v  in  the  turbulent 

m  m 

layer. 

Nondimens ional  variables  (wall  coordinates)  were  used  to  generalize  the 
analysis.  The  friction  velocity  is  defined  as 


*  /  /  x  1/2 

u  «  (x^/p) 

(2-2a) 

The  following  dimensionless  groups  were  formed: 

u+  *  u/u*. 

( 2-2b ) 

and 

y  ■  yu  /v  . 

(2-2c) 
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Equation  (2-1)  was  nond intensions l i zed  by  making  use  of  equations 
(2-2)  giving 


x/x  »  (1  +  e  /v)  du+/dy+  . 
w  m 


(2-3) 


Oeissler's  first-order  empirical  equation  for  eddy  diffusivity  of 
momentum  in  dimensionless  form  is 


c  /v  *  n2u+y+[l-exp(-n2u+y+)]  for  (0.0  <  y+  <  yt)  . 

ID  f 


(2-4) 


The  u+(y+)  relationship  for  the  sublayer  is  then  the  combination 
of  equations  (2-3)  and  (2-4),  which  gives  the  first-order  ordinary 
differential  equation 


du+/dy+  * 


x/x. 


1.0+n2u+y+[l.0-exp(-n2u+y+)] 


(sublayer)  . 


(2-5) 


This  relationship  holds  well  for  both  plain  tube  and  parallel  plate  channel 
+  2  ♦ 

flows  if  the  y^Ole^)  and  n  (ymt )  dependence  from  reference  15  is  used 

(figure  2-4).  A  polynomial  fit  of  these  data  was  made  so  that  the 
+  2  + 

y^(Re^)  and  n  (ymt)  dependence  shown  in  figure  2-4  could  be  included  in 
equation  (2-5).  See  tables  2-1  and  2-2. 

2 

The  index  n  in  equation  (2-5)  is  a  damping  factor  obtained  from  plain 

tube  experiments;  it  accounts  for  the  presence  of  the  wall  on  turbulence. 

2 

Including  this  n  (Re^)  relation  causes  the  sublayer  thickness  to  increase 
more  rapidly  with  decreasing  Re^  than  if  a  fixed  value  of  n2  was  used. 

This  implies  that  the  turbulent  viscosity  is  made  less  significant  compared 
with  the  molecular  viscosity  when  Re^  decreases  so  that  the  laminar 
u  -  y  relation  is  approached  (reference  15). 


Table  2-2.  Polynomial  Fit  of  Data  from  Table  2-1 

Refe  <  20,000  y/  »  -  0.691b  x  10"2  Reb  ♦  0.4103  x  10'6  (Refe)2  -  0.8668  x  10_U  (R«b>3  +  63.4614 

20,000  <  Re  <  50,000  y*  »  -  0.4339  x  10-3  Re  ♦  0.3822  x  10'8  (Re  )2  ♦  27.1373 
b  b  b 

Reb  >  50,000  yj  «  15.0 


y*  < 

205 

2 

n 

■  0.1408  x 

IQ'3 

♦ 

y.c 

-  0.2922  x 

10~* 

(y*  )2  -  0.0038 

■C 

205  <  y*. 

7  rat 

<  510 

2 

n 

>  0.1886  x 

io-4 

♦ 

y«t 

-  0.1709  x 

1 

O 

(y*)2  ♦  0.1813  X  io'u  (y*  )3  ♦  0.01003 

St  Bt 

y*  > 

510 

2 

n 

•  0.0154 

The  Prandtl  mixing  length  in  conjunction  with  Von  Karman'a  similarity 
hypothesis  yields  the  eddy  diffusivity  of  momentum  expression  for  the 
turbulent  layer  where 

t  •  k  (dY**~  2  \  '  cm  *  *2  (du/dy)*  and 
yd  u/dy  J 

k  a  0.36  (Von  Krfrmin's  constant), 

giving 


Substitution  of  equation  (2-6)  into  equation  (2-3)  gives  the  following 
u+(y+)  relation  for  the  turbulent  layer: 


. 2  + 
d  u 


+  x2 


dy 


+2 


(du*/dy*) 


(x/x  -du+/dy+)^2 


■  K 


(turbulent  layer) 


(2-7) 


The  value  0.36  for  Von  Karman's  constant  was  shown  to  yield  good 
results  for  the  limiting  cases  of  parallel  plate  channel  and  pipe  flow 
(reference  15). 

For  the  region  "i"  sublayer,  0.0  <  y^  <  ,  n  *  n^,  tw  * 

y+  *  y*  .  Equation  (2-5)  becomes 


du 


li 


t/t. 


.+  ,2+  +  ,  ,2++. 

dyu  ^Vliyli  l-«p(-ai“liyi1) 

where  u*^  »  0.0  at  y^  *  0.0  (from  boundary  condition  no.  2). 
turbulent  layer,  y  <  y^+<  y^,  equation  (2-7)  becomes 


(2-8a) 


In  the 


.2  + 

d  un 


dy 


+2 

2i 


+  >1/2 


where  uj£  (y^)  -  (y*,) 


(2-8b) 


Equations  (2-8)  were  integrated  to  obtain  the  velocity  profiles  for 
region  ”i”  using  a  fourth-order  Runge-Kutta  routine  as  described  in 
reference  18. 


Equation  (2-8a)  is  to  be  integrated  from  y^ 


0.0  to 


yli  "  yi/  ‘ 

The  step  slsea  H  -  yJ/ITER  and  H  ■  y*t/ITER2  were  used  for  the  sublayer  and 
turbulent  layer  calculations,  respectively.  A  parametric  study  of 
ITER2(Re^,b)  was  performed  to  minimize  computational  time;  from  this 
study,  figure  2-5  was  generated  and  a  curve  fit  determined.  The  expression 
for  ITER2  proved  to  yield  accurate  velocities  at  y  ■  y  when  compared  with 
velocities  resulting  from  very  large  values  of  ITER2  (see  '’Results"  later  in 
this  section). 
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Figure  2-5.  Suggested  Number  of  Steps  (in  Runge-Kutta  Routine)  for 
Turbulent  Region  Calculations 

At  the  end  of  the  integration  of  equation  (2-8a),  a  u^(y^)  profile 
would  have  been  generated  up  to  y^.  At  y|^  »  y^,  the  velocity  u^  and  the 
slope  of  the  velocity  profile  du^/dy^  oust  be  used  as  the  initial  values 
to  integrate  equation  (2-8b). 

A  similar  nondimens localising  procedure  was  used  to  obtain  differen¬ 
tiable  equations  necessary  to  produce  velocity  profiles  in  region  "o." 
Equation  (2-5)  becomes  (for  the  sublayer  0.0  <  y|Q  <  y*^) 


(2-9a) 


Uno  "to  'to  [1-“p(-”o  “to  yto)] 
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For  the  turbulent  layer  y+,  <  y+0  <  y*mt,  where  u+’  (yJJ  -  u[’  (y^)  , 
and  equation  (2-7)  become* 


.2  *' 

d  u2o 


dy 


+2 

2o 


( 2— 9b ) 


+ 1  + 

For  simplicity  u^Q  was  initially  set  equal  to  zero  at  y^Q  *  0.  After 

equations  (2-9)  were  integrated,  the  dimensionless  tube  velocity  was 

accounted  for  by 


*  u 


+  ' 
2o 


+  U 


+ 

tube 


The  Runge-Kutta  technique,  used  with  equations  (2-8),  was  similarly 
used  to  integrate  equations  (2-9)  for  region  "o." 


DERIVATION  OF  SHEAR  STRESS  RATIOS 

Equations  (2-8)  and  (2-9)  were  expressed  in  terms  of  the  shear  stress 

ratios  x/x.  and  x/x  .  These  ratios  vary  across  their  respective  boundary 
i  o 

layers.  A  force  balance  across  the  annular  gap  was  performed  to  find  an 
expression  for  these  .ratios.  This  may  be  done  only  for  the  fully  developed 
region  since  the  flow  is  not  accelerating.  The  following  derivation  was 
performed  to  find  how  these  shear  stress  ratios  vary  with  distance  from  the 
walls. 


The  continuity  equation  in  cylindrical  coordinates  is 


where  ux,  uf,  uQ  denote  the  partial  derivatives  of  the  axial  velocity  compo 
nent  with  respect  to  x,  r,  9,  etc.  The  quantities  v  and  w  represent  the 
radial  and  circumferential  velocity  components.  Applying  assumption  no.  4 
to  the  above  equation  gives 


CONTINUITY  ux  -  0  .  (2-10) 

The  linear  momentum  equations  in  cylindrical  coordinates  are 

r  MOMENTUM  J|*(0.0)v-p.-i|jtg.v  ^2v  - 

|l  +  (V  •  S)u  -  -  ^  |*  +  g  +  v?2u  . 

Applying  assumptions  1,3,  and  4,  and  neglecting  gravity  terma , 
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Integrating  equation  (2-llc)  twice, 

°  *  k  S  (r  *  Vv  *  °2)- 

Applying  boundary  condition  no.  1  across  the  gap,  r,  <  r  <  r  . 
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(2-12) 
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Now  applying  boundary  condition  no.  2  for  region  "i"  where  r  <  r 
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but  may  alao  be  found  from  boundary  condition  no.  3  for  region 
"o'*  where  r  >  - 
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To  simplify,  aaaume  Utub<  ■  0.0, 


C 


2 


mt  . 

—  lo,et( 


(2-13c) 


Equationa  (2-13b)  and  (2-13c)  were  set  equal,  thus  eliminating 
and  resulting  in  the  following  expression  for  r  : 
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To  obtain  an  expression  for  the  local  velocity  in  region  "o"  where 
r  >  r  equations  (2-13a)  and  (2-13c)  were  substituted  into  equation 

flit 

(2-12),  resulting  in 
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u  »  r^—  ^  fr^  -  r^  +  2r^  log  (r  /r)l  . 
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From  the  definition  of  average  velocity. 


(2-15) 
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Equation  (2-15)  was  substituted  into  equation  (2-16),  and  that  equation 
was  simplified  to 


S  A-*l  t  • 


(2-17) 


Then  the  local  velocity  in  terms  of  the  average  velocity  for  region  "o" 
was  obtained  by  combining  equations  (2-15)  and  (2-17),  eliminating  the 
dp/dx  term. 
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Using  the  definition  of  shear  stress,  x  ■  -  vrrj  ,  equation  (2-18)  becomes 
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The  ratio  of  local  to  tube  shear  stress  can  be  found  by  using 
equation  (2-19)  and  the  fact  that  x  ■  xq  at  r  »  rQ  : 
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(2-20a) 


Equation  (2-20a)  must  be  expressed  in  terms  of  the  dimensionless 

variables  y\  y*  ,  r\  so  that  they  may  be  inserted  into  equation  (2-9), 
o  rao  o 

which  can  then  be  integrated. 

Using  the  fact  that 
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and  the  relations  given  by  equations  (2-2),  where  tw  ■  x  ,  equation  (2-20a) 
is  transformed  to  the  dimensionless  form 
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A  similar  procedure  was  followed  to  obtain  the  dimensionless  expression 
for  T/tj. 

The  local  velocity  for  region  "i"  was  found  by  combining  equations 
(2-12),  (2-13a),  and  (2-13b),  giving 


4pv  %  ['!  -  ri  ♦  241°«.(ti/t)]  • 
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From  equations  (2-16)  and  (2-21),  the  average  velocity  is 
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Equations  (2-21)  and  (2-22)  were  combined  so  that  the  local  velocity  u 
could  be  expressed  in  terms  of  and  the  dp/dx  term  eliminated: 
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From  the  definition  of  shear  stress*  t  *  -  ydu/dr,  and  setting 
t  *  at  r  *  r^,  equation  (2-23)  yields 
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Using 
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(2-25b) 


and  equations  (2-2)  where  tw  *  for  region  "i,"  the  following  relationship 
was  obtained  for  region  "i"  in  terms  of  the  inner  variables; 
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The  ratio  to/t^  was  found  by  first  performing  a  force  balance  over  the 
entire  annular  region  and  then  on  region  ”i."  One  cannot  just  simply  use 
equations  (2-20d)  and  (2-25c)  to  find  since  they  were  found  by  consid¬ 

ering  regions  "i"  and  "o"  separately.  A  force  balance  across  the  annular  gap 
represented  by  the  area  enclosed  within  the  dashed  line  of  figure  2-6  yields 
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A  force  balance  across  region  "1"  represented  by  the  shaded  area  in 
figure  2-6  gives 


2Vi“  -  (rmt  -  ri)dp/dx 


(2-26b) 


The  dp/dx  term  was  eliminated  between  equations  (2-26 a)  and  (2-26b)(  and 
the  resulting  expression  for  was  reduced  to  the  following  by  making  use 

of  the  definitions  of  a  and  b  (see  "List  of  Symbols"): 
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(2-26c) 


From  equations  (2-25b),  (2-26c),  and  the  definitions  of  a  and  b,  the 
following  expression  for  y*^  in  terms  of  a,  b,  and  r*  was  found: 
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From  the  dimensionless  form  of  equations  (2-14)  and  (2-20c), 


y*  »  (1  -  a/b)r*  . 
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The  definitions  of  b,  r^,  r*,  and  equation  (2-26c)  were  used  to 
obtain  the  following  expression  for  r^  in  terms  of  r*: 


(2-28) 
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From  equation  (2-2a),  the  dimensionless  velocities  at  r  *  r  are 
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The  following  relationship  was  found  from  boundary  condition  no.  U  and 
equations  (2-30)  and  (2-26c): 
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DERIVATION  OF  REYNOLDS  NUMBER  FOR  CONTINUITY  CHECK 

To  ensure  constant  volume  flux*  the  bulk  velocity  or  Reynolds  number 
for  the  generated  velocity  profile  at  x  »  L*  must  match  its  value  at  the 
throat  of  the  annulus. 

From  the  definition  of  bulk  velocity  and  equations  (2-2), 
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The  quantities  U.  ,  and  U.  are  the  bulk  velocities  within  region  "i"  and 
region  "o,"  respectively  (figure  2-7).  The  quantities  u^(r  )  and  uQ(r  )  are 
determined  from  the  integration  of  equations  (2-8)  and  (2-9),  respectively. 

From  the  definition  of  Reynolds  number. 
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Part  1 


Part  2 


Equations  (2-25)  were  used  to  convert  part  1,  and  equations  (2-20)  were 
used  to  convert  part  2,  to  local  wall  variables  and  y*  .  Toe  following 
expression  resulted: 
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Equation  (2-34)  was  used  to  compute  the  bulk  velocity  in  order  to  per¬ 
form  the  continuity  check  and  ensure  constant  volume  flux.  Due  to  the 

smaller  size  of  U.  relative  to  Re.  ,  U.  was  used  to  check  continuity 
D  DO 

using  the  following  expression: 
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DERIVATION  OF  FRICTION  FACTOR  AND  DRAG  EQUATIONS 

The  goal  is  to  determine  the  drag  force,  which  acts  to  oppose  the  motion 
of  the  vehicle.  However,  as  an  intermediate  step,  friction  factors  are 
generally  used  as  a  means  of  reporting  both  analytical  and  experimental 
results. 

The  Fanning  friction  factor  f^  is  defined  as  follows: 
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Fp  =  drag  force  due  to  wall  shear. 


Ay  =  wetted  area, 
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This  friction  factor  may  be  expressed  nondimensionally  using  the  following 
relationships  in  conjunction  with  equation  (2-26c): 
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The  result  is 
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The  pressure  drop  was  calculated  by  first  finding  the  head  loss  h^  by 
using  the  Darcy-Weisbash  equation,  which  is  valid  for  laminar  or  turbulent 
duct  flows  of  any  cross  section: 


(2-38) 
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The  quantity  f^  is  the  Darcy  friction  factor,  which  can  be  found  from 
the  relation 


£d  *  4£f  • 


The  steady  flow  energy  equation  with  no  heat  transfer  or  shaft  work  and 
nondeformable  surfaces  for  incompressible  flow  through  a  constant  area  duct 
reduces  to 


(2-39) 


Equations  (2-38)  and  (2-39)  were  combined  to  give 
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The  total  drag  force  for  the  fully  developed  region  comprised  of  both  shear 
and  pressure  forces  is 


vSwWSS1 


where  L  is  the  entrance  length  and  Ly  is  the  vehicle  length.  The  quantities 
x.  and  tq  were  found  by  the  iteration  procedure  described  in  appendix  A. 

Once  the  shear  stresses  have  been  foundt  dp/dx  may  be  computed  from  equation 
(2-40)  and  the  relationship  between  the  Darcy  and  Fanning  friction  factors. 

The  shear  stress  and  pressure  drop  are  constant  in  the  fully  developed 
region;  therefore,  equation  (2-41)  reduces  to 

Dfd  ■  '  [2Vi  (Lv  -  L*)  -  ri  £  (Lv  -  L*)]  •  (Z-UZ) 


RESULTS 

A  computer  code  was  developed  following  the  iterative  scheme  described 
in  appendix  A.  A  parametric  study  was  performed  of  ITER2(Rev,b).  As  was 
previously  mentioned,  ITER2  was  the  parameter  that  controlled  the  step  size 
used  in  the  Runge-Kutta  routine.  For  each  Rev  and  b,  small  values  of 
ITER2  were  chosen,  then  increased  until  further  increase  in  ITER2  produced 
no  further  change  in 'the  computed  friction  factors  compared  to  a  very  large 
value  of  ITER2.  The  values  for  ITER2  that  resulted  in  +  3.0-percent 
difference  from  the  largest  ITER2  used  were  plotted  in  figure  2-5.  A  curve 
fit  performed  for  that  plot  is  also  displayed.  This  curve  fit  was  used  in 
the  computer  code  to  optimize  the  step  size  H  for  any  combination  of  b  and 

Re  .  Central  processing  unit  (CPU)  time  varied  from  10  minutes  to  4 

v  8 
hours.  The  longest  CPU  time  occurred  for  b  ■  2  at  Rev  »  10  . 

Figure  2-8  shows  the  convergence  to  a  uniform  velocity  profile  for  the 
case  b  ■  1.333  and  Ufc^  *  200.0  ft/s.  ITER#  ■  1  represents  the  first 
velocity  profile  computed  for  an  initial  guess  of  a  ■  .  This  is  a 


reasonable  starting  value  because  the  point  of  maximum  velocity  will  occur 

close  to  the  center  of  the  annulus.  Notice  that  for  ITER#  «  1,  there  exists 

a  large  discrepancy  between  the  values  of  U  (y  )  and  U  . (y  .). 

mo  mo  mi  mi 

This  discrepancy  is  a  physical  impossibility  and  served  as  a  measure  of  the 

error  of  the  first  guess  for  "a."  The  new  value  for  "a"  was  computed  by 

making  use  of  the  percentage  difference  between  U  .  and  U  to  steer  the 

mi  mo 

convergence  in  the  proper  direction.  Notice  at  ITER#  *  2,  r  .  was 

me 

increased,  which  effectively  pushed  the  region  "i"  profile  out  beyond  the 
region  "o"  profile  and  resulted  in  an  overcorrection.  However,  notice  that 
and  U  are  now  much  closer.  This  procedure  of  iterating  on  the 
quantity  a  resulted  in  rapid  convergence  to  a  continuous  velocity  profile  in 
just  five  iterations.  Remember  that,  although  a  continuous  velocity  profile 
has  been  found,  the  correct  physical  solution  is  not  guaranteed.  The 
profile  must  be  checked  for  continuity,  and  the  value  for  "a"  adjusted  in 
the  proper  sense  until  this  flow  condition  is  satisfied.  See  appendix  A. 


Figure  2-9  depicts  fc.  for  various  values  of  b  and  Re  .  The  data 

ta  v 

for  this  figure  are  in  appendix  B.  The  quantity  ff .  was  calculated  by 

ta-  2 

multiplying  f^  from  equation  (2-37)  by  the  factor  (b  -  l)/b  .  This 
factor  comes  from  the  continuity  relation,  equation  (A— 1 ) .  The  newly  found 
friction  factor  f^d  is  now  in  terms  of  Uv,  which  is  usually  known.  A 
considerable  amount  of  time  was  devoted  to  generating  this  figure  especially 
for  large  values  of  b  and  Rev,  which  required  a  large  number  of  iterations 
(see  figure  2-5).  Friction  factors  were  smallest  for  large  values  of  b  and 
Rev.  One  may  notice  from  figure  2-9  that  a  decrease  of  gap  width  by  a 
factor  of  2  results  in  increased  friction  factors  by  approximately  4000  to 
4500  times,  depending  upon  the  choice  for  Rey. 

The  following  relationships  were  determined  from  polynomial  regression 
of  observed  trends  throughout  the  Reynolds  number  and  radius  ratios  tested. 

The  value  for  the  quantity  a  was  found  to  increase  substantially  with  b 
and  only  slightly  with  Rev: 


a  -  Cl  [loglQ  (R«v)]  +  1.0 
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where 


Cl  *  0.001  (516.3b  -  516.737) 


for  1.01  <  b  <  1.143, 


C2  »  0.001  exp  (-  272.4b  +  607.9b  -  334.9) 


+  0.16X, 


Cl  *  0.001  (260.2b4  -  176.2b  -  59.04) 


for  1.143  <  b  <  2.0, 


C2  ■  0.001  exp  (-  3.407b  +  11.74b  -  4.804) 


+  1.07X. 


The  ratio  can  be  readily  computed  from  equations  (2-43)  and 


(2-26c)  for  given  value*  of  b  and  Uv: 


Rey  -  104, 
Rev  -  105, 
R«v  -  106, 
Rev  -  107, 
Rev  -  108, 


5mi  “  407-154b  “  243.765. 

6*  *  2933.95b  -  1755.34. 

6*.  »  27,267.1b  -  18,137.5. 
8*t  -  23,6729b  -  161,249. 

5*  -  2,072,380b  -  1,437,230. 


(2-44a) 

(2-44b) 

(2-44c) 

(2-44d) 

(2-44e) 


Theae  fits  may  come  in  handy  during  the  analysis  of  developing  flow  where  the 
upper  limit  6*^  is  needed. 


Figures  2-10  and  2-11  were  generated  from  shear  stress  data  obtained 
from  the  analysis.  The  shear  stresses  were  nondimensionalized  by  the 
following  equations: 


2-31 


LWD 


<V 

*  fd 

2 

fd  1/2  plT 

*T 


<To> 


fd 


p°; 

8. 


(2-45) 


Shear  stresses  on  the  vehicle  wall  may  be  readily  extracted  by  using  figure 
2-10  and  the  above  relations.  It  is  then  an  easy  task  to  compute  the  shear 
drag,  equation  (2-41).  Appendix  B  contains  the  values  used  to  generate 
figure  2-9. 

Figure  2-12  displays  plots  of  wall  variables  for  regions  "o"  and  "i"  for 
the  case 


Uv  »  200.0  ft/s, 
b  »  1.333, 
r^  *  6.0, 

v  *  1.786  x  10-^  ft2/s  (for  air  at  97.8°F), 

which  was  also  analyzed  in  reference  17.  Although  the  two-layer  model  was 
used  in  this  study,  the  two  plots  in  figure  2-12  are  representative  of  some 
three-layer  model  characteristics  exhibiting  a  sublayer  where  u+  ■  y+, 
the  buffer  layer,  and  then  the  turbulent  layer.  Figure  2-13  shows  how  the 
eddy  diffuaivity  of  momentum  varies  across  the  annulus.  This  figure  also 
shows  the  relative  importance  of  turbulent  to  molecular  viscosity  acrosE  the 
gap  of  the  core.  Both  figure  2-12  and  figure  2-13  show  close  agreement  with 
the  results  of  reference  17.  The  maximum  velocity  was  found  to  occur  at 
a  »  1.22138,  which  corresponds  to  66  percent  of  the  gap  width  measured  from 
the  vehicle  surface.  This  also  can  be  seen  in  figure  2-8  for  ITER#  *  5  and 
in  figure  2-13. 
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Figure  2-13.  Eddy  Diffuslvity  of  Momentum  Across  the  Annular  Gap 

CONCLUSIONS 

The  friction  factor  results  presented  in  figure  2-9  agree  with  the 
results  of  reference  19  within  4  to  13  percent,  and  with  the  results  of 
reference  S  within  approximately  12  percent.  The  fully  developed  friction 
factors  of  figure  2-9  were  compared  with  the  Moody  friction  factor  chart. 
This  comparison,  discussed  in  detail  in  appendix  C,  resulted  in  the  con¬ 
clusion  that  the  Moody  diagram  can  be  used  to  predict  the  pressure  drop  for 
fully  developed  turbulent  flow  providing  the  following  correlation  is  used: 
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where 


fM  i  0.001  +  2.8/(log10Reb) 


3.1 


As  seen  from  the  tabulated  results  in  appendix  B,  this  Moody  approximation 

is  accurate  to  within  4  percent  when  the  value  of  b  is  less  than  1.2  and 

has  errors  up  to  10  percent  for  b  between  1.2  and  2.0.  The  largest  error 

4  5 

occurs  at  lower  Reynolds  numbers,  10  -10  . 


3.  ANALYSIS  OF  DEVELOPING  FLOW 


EXPERIMENTAL  SURVEY 

The  region  of  developing  flow  is  the  part  of  the  flow  where  the  fluid 
core  accelerates,  and  therefore  the  velocity  profiles  are  functions  of  axial 
position.  The  acceleration  is  a  result  of  two  events. 

1.  The  fluid  at  a  station  just  ahead  of  the  vehicle  is  initially  at 
rest.  When  the  vehicle  approaches  this  station,  it  displaces  the  fluid, 
resulting  in  a  fluid  velocity  greater  than  its  initial  value  of  zero. 

2.  The  merging  of  boundary  layers  from  both  tube  and  vehicle  walls 
retards  the  flow  within  these  layers  and  causes  the  fluid  core  velocity  to 
increase  so  as  to  keep  the  volume  flux  constant  at  every  axial  location. 

At  the  annular  entrance,  the  boundary  layers  are  in  their  infancy.  The 
viscous  effects  have  not  had  sufficient  chance  to  propagate  and  are  confined 
to  extremely  thin  layers  next  to  the  walls.  The  velocity  of  the  fluid  core 
must  be  brought  to  the  respective  wall  velocity  in  a  very  short  distance. 
This  gives  rise  to  large  velocity  gradients  (du/dy)  at  the  walls.  These 
gradients  decrease  as  the  boundary  layers  continue  to  grow  downstream.  From 
the  relationship  xw  *  -p(du/dy),  it  becomes  clear  that  large  shearing 
stresses  are  present  when  large  velocity  gradients  exist,  especially  at  the 
annular  entrance.  It  is  important,  then,  to  determine  how  the  shear 
stresses  vary  with  axial  position  so  that  the  shear  drag  over  the  entire 
developing  portion  may  be  accurately  computed.  The  end  of  the  hydrodynamic 
entry  region  is  reached  when  the  upper  and  lower  boundary  layers  meet 
(figure  2-2). 

Unlike  fully  developed  annular  tube  flow,  little  analytical  and  experi¬ 
mental  information  is  available  on  turbulent  entry  flows  through  annuli, 
particularly  in  a  case  in  which  one  of  the  walls  is  moving. 

Experimental  investigations  of  turbulent  entry  flows  in  annuli  with 
stationary  boundaries  have  been  reported  in  references  20,  21,  and  22. 


In  reference  8,  the  investigators  measured  drag  forces  on  the  solid  core 
of  a  vertically  positioned  annulus.  The  core  was  suspended  from  a  calibrated 
spring  while  water  was  passed  downward  through  the  annular  gap.  Radius 
ratios  of  b  =  1.78  and  b  =  2.97  were  used  for  a  Reynolds  range  of  900  -  4.5 
x  10  ,  where  the  Reynolds  numbers  were  determined  from 


Re 


2(r0-4)gb 

r  v 
o 


The  length  of  the  annulus  was  varied  from  4  to  12  ft  so  that  entrance 
effects  could  be  studied.  Although  their  experimental  apparatus  was  crude, 
some  useful  facts  were  extracted.  The  increment  of  shear  friction  was  found 
to  become  persistently  larger  with  increased  Reynolds  numbers  as  measurements 
were  taken  closer  to  the  annular  entrance.  The  outer  wall  friction  factors 
were  computed  by  using  inner  wall  values  in  a  relationship  derived  from  a 
force  balance  of  the  annular  fluid.  This  relationship  may  be  used  only  for 
fully  developed  flow  since  a  force  imbalance  actually  exists  in  the  develop¬ 
ing  portion  of  the  flow.  The  outer  wall  factors  were  found  to  agree  within 
2  percent  with  the  data  in  a  previous  report  (reference  8)  for  the  range 
5000  <  Re  <  3.5  x  10^.  It  was  also  found  that  the  outer  wall  factors  were 

independent  of  the  radius  ratio.  In  addition,  the  r  ..(Re)  behavior,  which 

mt 

was  also  calculated  using  the  force  balance,  agreed  within  2  percent  with 

the  previous  data  for  Re  >  6000.  In  fact,  a  qualitative  assessment  of  the 

nature  of  the  f Iowl was  not  made  to  indicate  at  what  point  the  force  balance 

may  be  applied.  It  was  noted  only  that  beyond  250  equivalent  diameters  from 

the  throat,  the  data  were  little  affected  by  Reynolds  numbers.  Reference  8 

did  mention  that  the  use  of  the  force  balance  to  calculate  r  .  may  not  be 

DC 

appropriate  near  the  annular  entrance. 


In  reference  21,  experiments  were  conducted  on  the  entrance  region 
taking  static  pressure  measurements.  Reynolds  numbers  based  on  D.  and 
ranged  from  1.6  x  10  to  7  x  10  .  Two  radius  ratios  were  used, 
b  »  2.0  and  b  *  3.2.  Round  and  square  entrances  were  tested  (figure  3-1)  for 
each  of  these  ratios  to  study  the  effect  of  the  entrance  shape  on  developing 
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Figure  3-1.  Entrance  Configurations  Used  in  Reference  21 

flow.  Entrance  lengths  for  both  round  and  square  configurations  were  found 
to  be  20  to  25  equivalent  hydraulic  diameters,  respectively,  an  order  of 
magnitude  10  times  smaller  than  that  reported  by  reference  8.  The  entry 
length  for  turbulent  flow  in  circular  tubes  or  parallel  plate  channels  is 
usually  on  the  order  of  20  hydraulic  diameters.  The  investigators  in 
reference  21  found  that  the  square  configuration  caused  an  initial  separation 
period  in  which  the  pressure  drops  decreased  with  distance  from  the  throat. 

After  this  period,  pressure  drops  began  to  increase  for  approximately  three 
hydraulic  diameters  (figure  3-2,  a  and  b),  at  which  point  they  decreased 
again.  The  increase  in  the  pressure  drops  was  attributed  to  the  stream 
expansion  following  the  vena  contracta  downstream  of  the  sharp  entrance. 

For  the  rounded  entrances,  pressure  drops  were  found  to  decrease  up  to  1 

approximately  12  and  9  hydraulic  diameters  for  b  *  2.0  and  b  *  3.2, 
respectively,  then  increase  up  to  fully  developed  values  (figure  3-2c). 

The  inflection  points  for  both  of  these  ratios  were  attributed  to  transition 
from  laminar  to  turbulent  boundary  layers.  The  entry  length  for  the  rounded 
entrance  was  reduced  to  15  hydraulic  diameters,  and  a  monotonically 
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decreasing  pressure  gradient  resulted  by  tripping  the  boundary  layer 
turbulent.  The  rounded  entrances  were  found  to  generally  improve  flow 
stability.  The  results  shown  in  figure  3-2c  were  reported  to  be 
qualitatively  similar  to  those  in  reference  23  for  turbulent  flow  through 
rectangular  and  triangular  ducts  with  rounded  entrances. 

The  experiment  in  reference  23  seemed  to  be  carefully  done.  Reasonable 
entrance  lengths  comparable  to  plain  tube  and  parallel  plate  channel  flows, 
which  are  the  limiting  cases  of  annular  tube  flows,  were  obtained. 

The  sharp  entrance  condition  (figure  3-la)  is  analogous  to  the  case  in 
which  a  vehicle  is  just  exiting  a  tube.  The  smooth  entrance  case  (figure 
3-lb)  is  more  representative  of  the  condition  of  this  present  study,  except 
that  the  length  of  tube  preceding  the  core  may  be  insufficient.  The  data 
presented  in  reference  21  provide  good  insight  into  the  nature  of  developing 
annular  flows.  The  data  can  provide  a  baseline  for  comparison  with  an 
analysis  of  a  vehicle  just  before  it  exits  a  tube. 

In  reference  22,  experiments  were  conducted  similar  to  those  of 
reference  21.  Two  annuli  having  radius  ratios  of  b  *  1.88  and  b  *  2.91  were 
investigated.  Both  annuli  were  fitted  with  round  and  square  entrances 
(figure  3-3).  The  following  information  was  determined: 

1.  Mean-velocity  profiles 

2.  Outer  wall  static  pressure  gradients 

3.  Boundary  layer  thickness  parameters 

4.  Local  friction  factors. 
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Re^  ranged  from  7  x  10  to  1.6  x  10  .  Axial  development  of 
velocity  profiles  for  b  *  2.91  for  both  entrance  geometries  is  displayed  in 
figure  3-4.  The  velocity  profiles  for  the  square  entrance  are  skewed, 
especially  near  the  throat,  as  shown  in  figure  3-4a.  The  profiles  become 
less  skewed  and  the  maximums  shift  toward  the  center  of  the  gap  with 
increased  distance  from  the  annular  entrance.  The  profiles  for  the  round 
entrances  are  less  skewed  than  those  for  the  square  entrances,  looking  like 
slug  flow  near  the  entrance  (figure  3-4b).  The  shapes  of  the  profiles  were 
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Figure  3-3 .  Entrance  Configurations  Used  in  Reference  22 

found  to  be  independent  of  Reynolds  numbers  for  both  entrance  geometries. 
Entry  lengths  for  the  square  throats  with  b  »  1.88  and  b  »  2.91  were  found 
to  be  27.73  and  19.81  hydraulic  diameters,  respectively.  However,  entrance 
lengths  with  the  round  throats  were  unobtainable  because  of  flow 
instabilities  that  delayed  the  fully  developed  condition  beyond  the  test 
lengths.  This  contradicts  the  results  in  reference  21,  which  indicated  that 
the  entry  lengths  actually  decreased  from  25  to  20  equivalent  diameters  by 
changing  the  inlet  from  sharp  to  smooth.  Notice,  however,  that  the  inlet 
conditions  described  by  these  two  investigators  differ  in  that  for  reference 
22  the  core  preceded  the  tube  entrance.  Dissimilar  results  are  therefore 
not  unexpected. 

The  entrance  condition  investigated  in  reference  22  is  comparable  to 
that  of  a  vehicle  exiting  a  tube.  While  this  condition  is  not  exactly  that 
of  the  present  study,  the  results  of  reference  22  provide  useful  information 
on  developing  velocity  profiles.  The  location  of  transition  from  the 
laminar  to  the  turbulent  boundary  layer  was  not  measured  directly.  It  was 
estimated  using  shape  factors  where  a  value  of  1.4  (from  flat  plate  theory) 
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was  used  to  indicate  transition  from  the  laminar  to  the  turbulent  boundary 
layer.  The  calculated  shape  factors  never  reached  a  value  of  1.4. 
Furthermore,  the  shape  factor  that  indicates  transition  may  be  something 
other  than  1.4  for  annular  tube  flow.  Thus,  it  is  not  certain  where 
transition  occurred,  if  at  all.  Nevertheless,  the  investigators  contend  that 
transition  occurred  near  the  throat  of  both  annuli.  Pressure  drop  data  for 
both  entrance  geometries  showed  similar  behavior  to  the  data  in  reference  21. 

More  recently,  the  experiments  reported  in  references  6,  24,  25,  26,  and 
27  were  aimed  at  obtaining  correlations  between  drag  force  on  a  vehicle, 
blockage  ratio,  and  Rev-  Davidson,  in  reference  6,  provided  a  brief  review 
of  these  experiments  together  with  a  compilation  of  data  onto  a  single  curve 
(figure  3-5).  The  experimenters  chose  0.5  x  10J'for  Rey.  Davidson  also 
presented  pressure  signatures  within  the  annulus  from  his  experiments. 
was  predicted  from  our  present  theory  for  various  nose  coefficients  C  , 
and  plotted  on  this  figure.  A  pictorial  description  of  nose  shapes  and  their 
associated  nose  coefficients  are  displayed  on  page  3-12  of  reference  30. 

ANALYTICAL  SURVEY 

Analytical  attempts  to  model  the  entry  region  were  reported  in 
references  4,  5,  and  22.  Only  round  entrances  may  be  considered  for 
conventional  boundary  layer  analysis  since  separation  occurs  in  flows  with 
square  entrances. 

In  reference  22,  it  was  assumed  that  radial  pressure  variations  and 

fluctuating  velocity  components  were  negligible.  The  momentum  integral 

equation  was  applied  to  both  inner  and  outer  boundary  layers  separately.  The 

fluid  core  outside  the  boundary  layers  was  assumed  irrotational.  The 

experimental  results  showed  that  the  irrotational  flow  field  existed  from  the 

annular  entrance  up  to  10-14  hydraulic  diameters  from  the  throat  for  the 
4  4 

Re^  range  7  x  10  to  1.6  x  10  .  Trends  in  local  friction  factors  were 
reportedly  similar  to  the  data  of  reference  28  in  which  friction  factors  were 
generated  by  using  pressure  drop  data  at  x/D^  >  3.81  and  x/Dh  >  5.33  for 
b  »  2.91  and  b  ■  1.88,  respectively,  where  the  momentum  flux  was  assumed 
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Figure  3-5.  Comparisons  of  Experimental  Drag  Coefficients  of  Principal 

Investigators  (from  Reference  6) 


negligible  as  observed  from  experimental  data.  The  momentum-integral 
equation  was  reduced  to  the  following  in  terms  of  friction  factors: 
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It  seems  odd  that  the  momentum  flux  becomes  negligible  so  soon  before 

fully  developed  flow,  which  occurred  beyond  x/D,  »  24.95  and  x/D,  *  34.93 

h  h 

for  b  »  2.91  and  1.88,  respectively. 

Wilson  and  Medwell  (reference  4)  also  used  the  momentum  integral  method 

for  radius  ratios  b  *  1.25  to  b  a  5.0  over  a  range  of  Re.  from  104  to 

5  “ 

3  x  10  .  Reichart's  expression  for  eddy  diffusivity  as  suggested  in 

reference  29  was  modified  by  multiplying  it  by  the  Van  Driest  damping  factor. 

This  modified  expression  produced  velocity  profiles  that  compared  well  with 

the  experimental  data  of  reference  22  (figure  3-6).  The  hydrodynamic  entry 

length  for  all  the  radius  ratios  and  Reynolds  numbers  considered  by  Wilson 

and  Medwell  was  found  to  be  approximately  10  equivalent  diameters,  which 

strengthens  the  idea  that  hydrodynamic  development  is  weakly  dependent  on 

radius  ratios  for  the  range  tested.  It  appears  that  Wilson  and  Medwell  have 

successfully  predicted  developing  profiles  in  accordance  with  the 

experimental  data  of  Okiishi  and  Servoy  (reference  22). 

Like  Wilson  and  Medwall,  Sud  and  Chaddock  (reference  5)  used  the 
momentum  integral  technique.  The  eddy  viscosity,  however,  was  modeled  using 
a  two-layer  approach  suggested  by  Quarmby  in  reference  2.  Deissler's 
equation  for  eddy  diffusivity,  equation  (2-5),  was  used  for  sublayer  calcu¬ 
lations,  and  Von  Kdrada’t  similarity  hypothesis,  equation  (2-7),  was  used  to 
obtain  the  eddy  diffusivity  expression  for  turbulent  layer  calculations.  Sud 
and  Chaddock  made  provisions  for  a  moving  boundary,  which  previously  had  not 
been  done.  Only  one  radius  ratio,  b  »  1.33,  and  two  tunnel  pressures,  0.1 
and  0.2  atmospheres,  were  considered  at  Re.  »  1.02  x  106  and  Re.  * 

2.55  x  10°,  respectively.  Sud  and  Chaddock  claim  that  their  entry  analysis 
for  large  x/D^  predicted  velocity  profiles,  wall  shear  stresses,  and  fric- 


tion  factors  that  compared  favorably  with  the  corresponding  values  predicted 
by  their  fully  developed  flow  analysis.  But  closer  examination  of  the  plot¬ 
ted  profiles  (figure  3—1 lb )  reveals  that  the  volume  flux  was  not  held 
constant  throughout  the  entry  region,  rendering  their  analysis 
questionable.  The  total  drag  coefficient,  however,  was  found  to  be  in  good 
agreement  with  the  compilation  of  reference  6  (figure  3-5).  For  the  generic 
case  of  no  moving  boundaries,  Sud  and  Chaddock  also  claimed  good  agreement 
with  the  data  of  Quarmby  in  reference  3. 

Self-consistency  and  reasonable  agreement  with  existing  data  convinced 

the  author  of  the  present  study  to  take  an  approach  similar  to  that  of  Sud 

and  Chaddock.  However,  more  care  was  taken  throughout  the  analysis  of 

developing  flow  in  keeping  the  volume  flux  constant.  The  upper  limit  of 

+  2 

integration  for  sublayer  calculations  y^  and  the  damping  factor  n  were 
not  held  fixed  as  suggested  by  Sud  and  Chaddock,  but  were  allowed  to  vary 
with  Re^  and  u+  as  suggested  by  Quarmby  in  reference  15  and  as  discussed 
under  "Assumptions  and  Boundary  Conditions"  in  section  2.  The  approach  used 
in  this  study  is  summarized  as  follows: 

1.  The  eddy  diffusivities  for  the  sublayer  and  turbulent  layer  were 
modeled  by  Deissler's  equation  and  Von  Harmin' s  similarity  hypothesis, 
respectively. 

2.  Velocity  profiles  were  generated  by  integrating  the  governing 
differential  equations  using  the  Runge-Kutta  technique.  The  profiles  were 
modified  iteratively  until  continuity  requirements  were  satisfied. 

3.  The  functional  dependence  of  x^  and  Ap  with  x  were  found  by 
implementing  the  momentusi-integral  technique  coupled  with  the  conservation 
of  mass  for  boundary  layers  "o"  and  "i"  and  for  the  whole  annular  gap.  The 
three  resulting  equations  were  solved  simultaneously  (as  discussed  later  in 
this  section)  for  the  unknowns  Ap  and  Ax. 

4.  Having  determined  tw(x)  from  step  3,  the  shear  drag  was 
calculated  by  numerically  integrating  the  shear  stresses  over  the  entry 
length  and  multiplying  the  result  by  that  part  of  the  vehicle's  surface  area 
located  within  the  developing  region. 

5.  From  step  3,  the  Ap(x)  relationship  was  determined.  These  pressure 
drops  were  numerically  integrated  over  the  entry  length  to  obtain  the  total 


pressure  drop  across  that  region.  The  profile  drag  was  then  found  by 
multiplying  the  total  pressure  drop  by  the  vehicle's  projected  area. 


ASSUMPTIONS  AND  BOUNDARY  CONDITIONS 

The  same  boundary  conditions  as  listed  in  section  2  apply.  The 
following  additional  assumption  was  made  for  the  entry  analysis. 

The  velocity  profiles  u(y)  were  determined  from  equations  (2-8)  and 
(2-9)  before  integrating  equations  (3-11),  (3-12),  (3-16),  (3-23),  and 
(3-24).  Unlike  the  analysis  for  the  fully  developed  flow,  the  shear  stress 
ratios  t/t^  and  t/to  were  not  determined  from  equations  (2-25c)  and  (2-20d) 
that  resulted  from  a  force  balance.  The  accelerating  fluid  in  the  develop¬ 
ing  flow  region  suggests  that  a  force  imbalance  exists.  These  ratios  were 
linearized  from  1  at  the  walls  to  0  at  the  edge  of  the  boundary  layers  as 
suggested  by  Wilson  and  Medwell  in  reference  4. 


LOCATION  OF  VELOCITY  PROFILES 


The  linear  momentum  equation,  coupled  with  the  conservation  of  mass 
principle,  was  used  to  determine  the  axial  position  of  each  velocity  profile 
generated.  The  upper  and  lower  boundary  layers  were  considered  separately. 
Furthermore,  two  control  volumes  were  used  for  each  boundary  layer.  The 
control  volume  used  to  locate  the  first  profile  within  the  boundary  layer  of 
region  "i"  is  shown  in  figure  3-7a.  Subsequent  profiles  in  region  "i"  were 
located  using  the  control  volume  displayed  in  figure  3-9. 


Consider  a  control  volume  for  region  "i"  boundary  layer  as  shown  in 
figure  3-7a.  For  a  steady,  fixed  control  volume,  the  conservation  of  mass  is 


n)dA  ■  0  , 


(3-1) 


where  V  is  composed  of  u,  v,  w,  the  axial,  radial,  and  circumferential 
velocities.  Equation  (3-1)  becomes 


p(V  •  n)dA  +  f  p(V  •  n)dA  +  f p(V  •  n)dA 


■p(V  •  n)dA  =  0  , 


(3-2) 


where  the  integrals  over  1,  2,  3,  and  4  represent  the  four  boundaries  of  the 
control  surface.  The  results  of  the  dot  products  for  this  control  volume  in 
the  above  equation  are  summarized  in  table  3-1. 

Table  3-1.  Resultant  Dot  Product  on  Control  Volume  Surfaces 


Section  No. 

Axial  Station 

Boundary 

Layer  Height 

Resultant 

Dot  Product 

1 

X  a  0 

a  h 

V 

•  n  a  -U. 

in 

2 

0  <  x  <  xl 

h  <  $i  <  h  +  $* 

V 

•  ii  a  0 

3 

X  a  xl 

5.  =  h  +  S*. 

l  l 

V 

•  n  a  u 

4 

0  <  x  <  xl 

S.  a  0 

1 

V 

•  n  a  0 

Applying  the  information  from  table  3-1  to  equation  (3-2), 
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The  steady  linear  momentum  equation  is 
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Applying  the  information  from  table  3-1  to  the  above  equation  yields 
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Insert  equation  (3-3)  into  equation  (3-4): 
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The  only  forces  on  the  control  volume  are  those  due  to  shear  and 
pressure  (see  figure  3-7b): 
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To  find  A  (the  surface  area  of  the  cone  created  by  the  streamline), 
s 

the  streamline  depicted  in  figure  3-7a  was  assumed  to  grow  linearly 
with  x  (see  figure  3-8): 


X1 

J  (H  +  x  tan  (3)dx  =  2irx^  £r^  +  6^  +  ( / 2 ) J  , 


A  =  2ir 
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and 


(A  +  AA) .  =  ir(rj .  -  r?)  . 
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can  be  determined  from  equation  (3-3): 
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which  simplifies  to 
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Figure  3-8.  Assumed  Linear  Streamline  Growth 
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Equation  (3-9)  is  a  quadratic;  then 


2(ri+5i)  + 


The  function  f(r)  =  (u/U.  -  l)r  is  bounded  between 

in 


-r^  <  f(r)  <  0  . 


Then  the  bounds  for 
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f (r)dr 


~(ri  +  &i)&.  <  M  <  0 


For  the  lower  limit  of  M  *  -(r^  +  S^)» 
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The  only  reasonable  solution  is  obtained  when  the  negative  root  is  taken 


s*  •  ri  *  h  -  <ri  ♦  h*  y 


For  the  upper  limit  of  H  ■  0, 


Si  *  ^ri  +  5i^  -  ^ri  +  5i^ 


Once  again,  the  only  reasonable  solution  is  obtained  when  the  negative 
root  is  taken,  giving 


S*  =  0. 
1 


Therefore,  for  all  computations,  the  negative  root  should  be  used, 

* 

resulting  in  the  following  expression  for  S.: 
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and  incorporating  equations  (3-5)  and  (3-6)  gives 
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Similarly,  for  region  "o"  the  momentum  balance  gives 

Region  "o" 


-  (A  +  AA)  Ap  -  irr  (x,  +  t.  )  x,  = 
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The  unknowns  in  equations  (3-11)  and  (3-12)  are  Ap,  x, ,  (t .  ).  and  (x,  )  . 

1  in  i  in  o 

Initially,  it  was  assumed  that  (x.  ).  *  (x, ).  and  (x.  )  *  (x, )  .  After 

m  i  it  mo  l  o 

several  values  for  x  (x)  were  found,  then  (x.  ).  and  (x.  )  were  found  by 

w  in  i  mo 

extrapolating  x^(x)  out  to  x  *  0.0.  Therefore,  the  only  two  unknowns 
present  in  equations  (3-11)  and  (3-12)  are  Ap  and  x^.  However,  since 
equations  (3-11)  and  (3-12)  do  not  sufficiently  describe  the  whole  flow, 
another  relationship  is  needed. 

Consider  a  control  volume  that  incorporates  the  whole  annular  gap  as 
displayed  in  figure  3-9.  Summing  the  forces  on  the  control  volume. 
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Figure  3-9.  Control  Volume 
of  Annular  Gap 


Applying  the  continuity  equation  across  surfaces  and  gives 
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The  momentum  equation  across  sections  and  gives 
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Combining  equations  (3-13),  (3-14),  and  (3-15)  and  simplifying  gives  the 
following  relationship: 
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The  two  unknowns  Ap  and  are  found  for  region  "i"  by  solving 
equations  (3-11)  and  (3-16)  simultaneously.  For  region  "o,"  equations  (3-12) 
and  (3-16)  are  solved  simultaneously.  Naturally,  the  values  of  Ap  and  x^ 
for  region  "o"  will  not  initially  agree  with  their  counterparts  from  region 
"i.”  The  iterative  process  described  in  appendix  D  must  be  performed  until 
the  region  "i"  and  region  "o"  predictions  of  x^  match.  As  a  consequence, 
the  Ap  from  both  regions  will  also  match. 

To  evaluate  the  location  of  the  velocity  profiles  within  region  "i" 
down  stream  of  the  first  profile,  a  momentum  balance  was  performed  between 
two  successive  profiles  as  depicted  in  figure  3-10.  The  linear  momentum 
equation  for  this  control  volume  is 


IFx  »  - pn" Ak  -  UgavAm  +  J" p(u  +  Au)2d(A  +  AA)  , 


where  U^av  is  the  average  core  velocity  between  x  *  xn  and  x  *  xn+^- 
Using  the  fact  that 
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equation  (3-17)  becomes 
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b.  PROFILE  2 


Figure  3-10.  Control  Volume  Used  to  Locate  Successive  Velocity  Profiles 

Within  Region  "i"  Boundary  Layer 


From  continuity, 
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The  sum  of  the  forces  on  the  control  volume  (see  figure  3-10b)  is 
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The  boundary  Layer  was  assumed  Co  grow  linearly  with  x;  thus,  the 
following  expression  can  be  written: 
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where 
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Inserting  equations  (3-20b)  through  (3-20e)  and  equation  (3-21)  into 
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equation  (3-20a)  and  neglecting  in  terms,  gives 
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Combining  equations  (3-18),  (3-19),  and  (3-22),  and  letting 
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the  following  expression  is  arrived  at  for  region  "i": 
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Then  dp  and  dx  can  be  found  for  region  "i"  by  solving  equations  (3-23)  and 
(3-16)  simultaneously.  Similarly,  dp  and  dx  for  region  "o"  can  be  found  by 
solving  equations  (3-24)  and  (3-16)  simultaneously.  There  is  no  guarantee 
that  the  prediction  of  dx  from  region  "i"  will  match  that  from  region 
"o."  The  quantity  S*  was  iteratively  changed  until  the  predicted  value 
of  dx  from  region  "i"  matched  that  from  region  "o."  Once  the  dx  predictions 
from  both  region  ”i"  and  region  "o"  equations  are  in  close  agreement,  the 
dp  prediction  from  both  regions  should,  as  a  consequence,  also  be  in 
close  agreement.  For  details  of  the  iterative  procedure,  see  appendix  D. 

The  irrotational  assumption  resulting  in  the  expression 

dp/dx  =  -  pUs(dU&/dx)  (3-25) 

was  not  used  within  the  momentum  balance.  Once  dp  and  x^  or  dx  were  found 
from  the  momentum  balance,  they  were  compared  with  dp/dx  from  the  irrota¬ 
tional  assumption.  This  was  done  to  determine  if  the  momentum  balance  was 
performed  correctly. 

SCHEME  FOR  DETERMINING  VEHICLE  DRAG  IN  DEVELOPING  FLOW 

The  scheme  shown  in  appendix  D  was  used  to  obtain  the  local  pressures 
and  shear  stresses  in  the  annulus  during  the  developing  portion  of  the  flow 
up  to  the  point  where  the  boundary  layers  meet  (fully  developed).  The 
scheme  provided  the  following: 

1.  Convergence  on  matching  of  and  UgQ. 

2.  Constant  volume  flux  throughout  the  annulus. 

3.  Correct  axial  position  of  the  developing  velocity  profiles  and 
boundary  layer  heights. 

The  general  procedure  for  obtaining  these  results  is  as  follows: 

1.  Equations  (2-8)  and  (2-9)  were  used  to  generate  the  velocity  profiles 
for  each  of  the  chosen  dimensionless  boundary  layer  heights  6?. 


2.  The  following  procedure  was  implemented  so  that  Ap(x)  and  t^tx) 
could  be  determined. 

a.  To  locate  first  velocity  profile: 

(1)  Solve  equations  (3-11)  and  (3-16)  simultaneously  for  x^  and 
Ap. 

(2)  Solve  equations  (3-12)  and  (3-16)  simultaneously  for  x^  and 
Ap. 

(3)  Change  5*  and  follow  iteration  scheme  of  appendix  D. 

(4)  Repeat  steps  (1)  through  (3)  until  x^  values  predicted  from 
steps  (1)  and  (2)  match. 

b.  To  locate  successive  velocity  profiles: 

(1)  Solve  equations  (3-23)  and  (3-16)  simultaneously  for  Ax  and  Ap 

(2)  Solve  equations  (3-24)  and  (3-16)  simultaneously  for  Ax  and  Ap 

(3)  Change  6*  and  follow  iteration  procedure  of  appendix  D. 

(4)  Repeat  steps  (1)  through  (3)  until  Ax  values  predicted  from- 
steps  (1)  and  (2)  match. 

3.  The  local  pressure  drops  and  shear  stresses  on  the  vehicle  wall  at 
each  axial  station  have  been  determined  from  step  2. 

The  fluid  dynamic  drag  imparted  on  the  portion  of  the  vehicle  subjected 
to  developing  flow  is 


The  total  drag  on  the  vehicle  due  to  concentric  annular  flow  (developing 
and  fully  developed)  is  obtained  by  summing  equations  (2-42)  and  (3-26). 

The  drag  contributions  due  to  the  nose  and  base  of  the  vehicle  have  not 
been  incorporated  in  the  analysis.  They  can  be  included  to  find  the  overall 
drag  as  discussed  in  appendix  H. 


RESULTS 


Fifteen  velocity  profiles  were  typically  generated  within  the  developing 

region.  The  CPU  time  was  seen  to  vary  from  10  to  40  minutes.  Generally, 

the  CPU  time  increased  with  increasing  values  of  Rev  and  b  for  0.01  -  b  5  2.0 

and  10^  <  Re  <  10^. 

-  v  - 

A  test  case  was  run  so  that  the  present  analysis  could  be  compared  with 
that  of  Sud  (reference  17).  Figures  3-lla  and  3—1 lb  show  this  comparison 
for  the  case  of: 

Utube  =  200  ft's 

r^  =  6  ft 

r0  *  8  ft 

p  *0.1  atm 

Reb  *  1.02  x  106 

T  *  98°F. 

Profiles  l  and  2  of  figure  3-lla  are  the  developing  profiles  at 

x/D,  *  1.53  and  x/D,  *  18.95,  respectively.  Profile  2  is  the  nearly 
h  u 

fully  developed  velocity  profile  as  predicted  by  the  developing  analysis. 
Profile  3  was  computed  from  the  fully  developed  analysis.  Note  that  profiles 
2  and  3  are  very  similar,  even  though  they  were  computed  independently,  with 
continuity  being  the  only  common  link.  Closer  examination  of  these  profiles 
reveals  that  the  area  under  each  curve  is  the  same,  which  indicates  that  the 
volume  flux  remained  constant  throughout  the  analysis.  The  near  equality 
between  profiles  2  -and  3  gives  confidence  that  the  current  model  is  self- 
consistent  and  correct  within  the  limits  of  integral  theory.  The  profiles 
of  figure  3-lla  and  those  obtained  by  Sud,  figure  3-llb,  are  similar  with 
respect  to  profile  shape,  magnitude,  and  position  of  maximum  velocity. 

Careful  examination  of  the  profiles  labeled  "x/D^  ■  11.^7"  and  "Fully 
Developed"  in  figure  3-llb  suggests  that  continuity  was  violated. 

Figure  3-12  shows  that  the  present  model  compares  well  with  the  experi¬ 
mental  results  of  Davidson  (reference  6).  The  pressure  drop  over  Ln  could 
be  computed  by  using  a  nose  coefficient  from  Hoerner  in  conjunction  with  the 
Bernoulli  equation  (see  appendix  H).  The  exact  coefficient  for  an  ogive 
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Figure  3-12.  Comparison  of  Experime ital  Results  from 
Reference  6  and  Present  Theory 
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nose  was  not  known;  therefore,  the  Ap  across  Ln  was  set  equal  to  the 
corresponding  value  obtained  by  Davidson.  The  present  analysis  takes  over 
from  that  point  on.  The  Ap  across  length  2  agrees  within  2  percent 
with  Davidson's  data  (see  figure  3-12). 

The  following  example  was  run  because  it  represents  a  typical  N’avy 
problem  of  a  torpedo  moving  within  a  tube. 

Utube  =  ft/s, 
r^  =  10 . 5  in. , 

r0  »  12.0  in. , 

v  =  0.12583  x  10“^  ft^/s  (sea  water  at  60°F). 

The  results  of  the  developing  flow  model  are  displayed  in  figures  3-13 
through  3-18.  Figure  3-13  shows  3  of  15  computed  velocity  profiles  as  they 
develop  in  the  annulus.  These  are  at  approximately  10  percent,  50  percent, 
and  95  percent  of  boundary  layer  growth.  Appendix  G  (table  G-l)  summarizes 
the  results  of  the  entry  computation  shown  in  this  figure.  Notice  that  the 
developing  flow  model  at  large  x/D^  agrees  well  with  the  fully  developed 
flow  model.  The  largest  discrepancy  (3.1  percent)  exists  for  the  Sq  computa¬ 
tion.  A  comparison  of  core  velocities  reveals  that  they  agree  within  1.3 
percent  for  this  case.  The  volume  fluxes  within  the  three  regions  (lower 
boundary  layer,  core,  and  upper  boundary  layer)  for  these  three  profiles  are 
listed  in  table  G-2.  It  is  interesting  to  note  that  throughout  development 
the  volume  flux  for  the  region  "i"  boundary  layer  is  always  greater  than 
that  for  region  "<JV"  ’  This  trend  was  consistent  throughout  the  range  of  b 
and  Re  tested.  The  total  volume  fluxes  Q  for  each  profile  displayed 
in  figure  3-13  is  within  0.8  percent  of  the  volume  flux  entering  the  annulus 
Q^.  This  indicates  that  continuity  has  been  fulfilled  throughout  the 
developing  flow  analysis. 

Figure  3-14  shows  shear  stress  development  on  both  the  tube  and  vehicle 
walls.  Notice  that  maximum  shear  stress  on  both  walls  occurs  at  the  throat 
of  the  annulus.  The  maximum  developing  shear  stress  on  the  vehicle  is 
approximately  7.7  times  greater  than  the  value  for  fully  developed  flow.  By 
approximately  4  hydraulic  diameters  (15  percent  of  entry  length)  the  shear 
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stresses  decay  to  within  10  percent  of  the  values  for  fully  developed  flow. 
The  rapid  decay  of  these  developing  shear  stresses  translates  to  a  small 
increase  in  shear  drag  when  compared  with  fully  developed  values. 

Figure  3-15  displays  the  local  pressure  drops  during  profile 
development.  The  pressure  drop  at  the  annular  throat  is  roughly  6.3  times 
greater  than  the  value  for  fully  developed  flow  and  decays  to  within  10 
percent  of  such  values  by  approximately  11  hydraulic  diameters  (41  percent 
of  entry  length).  Of  the  two  (shear  and  form)  drag  contributors  in  the 
entry  region,  the  form  drag  seems  to  be  of  greater  significance  to  the 
total.  To  quantify  these  two  drag  contributors,  the  drag  force  due  to  each 
(calculated  in  appendix  E)  is  presented  below: 


31.9  percent 


68.2  percent 


(PDRG)g  =  -14,788  lbf  [27.5  percent  of  (DRG)^] . 


(SHRDRG)g  =  -2354  lbf  [4.4  percent  of  (DRG)  ] . 


(SHRDRG)fd  *  -5424  lbf  [10.1  percent  of  (DRG)T1 . 


(PDRG)fd  =  -31,252  lbf  [58.1  percent  of  (DRG)^] . 


(DRG)T  =  -53,818  lbf. 


As  suspected,  (PDRG)fi  was  substantially  larger  than  (SHRDRG)^.  The 
largest  drag  contributor  is  (PDRG)^,  with  the  fully  developed  flow  region 
contributing  68.2  percent  to  the  total  drag.  Of  course,  the  amount  that  all 
of  these  contribute  to  drag  would  vary,  depending  on  Lv,  b,  and  Uy. 

Figure  3-16  shows  that  use  of  the  irrotational  assumption  within  the 
momentum  balance  can  be  an  inaccurate  method  for  computing  dp/dx.  The 
irrotational  assumption  does  rather  well  between  the  second  computed  point 
and  at  four  hydraulic  diameters.  Beyond  four  hydraulic  diameters,  dp/dx 
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HYDRAULIC  DIAMETERS  DOUN  THE  ANNULUS 

Figure  3-15.  Pressure  Signature  in  the  Annulus  for  Developing  Flow 
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Figure  3-16.  Pressure  Signature  of  Present  Method  Compared  with  That 

from  Irrotational  Assumption 


dips  below  the  fully  developed  value  long  before  the  flow  becomes  fully 
developed  (30  percent  of  entry  length). 

Figure  3-17  displays  the  friction  factor  behavior  during  flow  develop¬ 
ment.  It  takes  approximately  11  hydraulic  diameters  for  the  friction  factor 
to  decay  to  within  10  percent  of  values  for  fully  developed  flow.  The 
largest  factor  occurs  at  the  annular  throat  and  is  approximately  six  times 
greater  than  the  fully  developed  value. 

Figure  3-18  displays  the  velocity  ratio  Ug/Ug^  as  a  function  of 
hydraulic  diameters  down  the  annulus.  It  is  important  to  note  that  the  core 
velocities  are  increasing  monotonically .  This  behavior  must  occur  and 
acts  as  an  indicator  of  the  numerical  stability  of  the  analysis.  Notice 
from  this  figure  that  the  velocity  at  the  annular  throat  is  approximately  9U 
percent  of  the  fully  developed  value.  If  the  monotonicity  of  the  core 
velocities  is  not  achieved,  the  dUc/dx  may  fluctuate  wildly  and  cause 
considerable  error  in  calculations  for  the  local  pressure  drops.  This  can 
be  easily  seen  from  the  irrotational  assumption 


The  curves  of  figure  3-19  were  generated  by  computing  dp/dx  at  a  number 
of  axial  locations  (same  as  figure  3-15)  for  each  R«v  and  b.  The  inte¬ 
grated  value  was  their  determined  from 

* 

T 

b  /  <s>- 

0 

and  nondimensionalized  as  shown  in  the  upper  right  corner  of  figure  3-19. 

The  curves  of  figures  3-20  and  3-21  were  generated  in  a  manner  similar 

to  that  used  for  figure  3-19;  however,  the  values  for  ftjJ  andft  J  were 

'e  v  0Je 

nondimens ionali zed  as  shown  in  figures  3-20  and  3-21,  respectively. 


Figures  3-19  and  3-20  may  be  used  to  determine  the  drag  only  for  the 
developing  portion  of  the  flow  in  the  annulus.  However,  to  actually  calcu¬ 
late  the  drag  (shear  and  form)  in  this  region,  entry  length  L*  is  needed. 
Figure  3-22  pro%rides  this  information.  Appendix  E  shows  how  one  would 
compute  the  drag  force  on  a  vehicle  given  Uv,  Lv,  r^,  and  b  using  only 
figures  2-9,  2-10,  3-19,  3-20,  and  3-22.  Data  tables  used  to  generate  curves 
for  figures  3-19,  3-20,  3-21,  and  3-22  may  be  found  in  appendix  F. 

4 

Results  for  Re  =  10  could  not  be  obtained  because  of  the 
v 

following  reasons: 


4.  CONCLUSIONS  AND  RECOMMENDATIONS 


Several  conclusions  may  be  drawn  from  this  study: 

L.  The  analysis  presented  herein  has  been  proven  to  be  self-consistent . 
The  developing  flow  analysis  for  large  x/Dh  predicts  nearly  identical  fully 
developed  values  for  velocity  profiles,  x^,  xQ ,  S^,  and  Sq.  See  appendix  G. 

2.  The  friction  factors  for  fully  developed  flow  are  in  good  agreement 
with  the  results  of  references  5  and  19. 

3.  Assuming  the  annulus  to  be  entirely  fully  developed  would  result 
in  lower  drag  estimates  for  L^/d  >  12.  However,  for  blunt  vehicles  with 
Ly/d  <  12,  the  flow  will  be  largely  in  a  developing  state.  Thus,  a  fully 
developed  approximation  would  yield  extremely  low  drag  predictions  if  applied 
toward  the  whole  flow. 

4.  The  Moody  chart  method  (excluding  entrance  and  exit  losses)  can  lead 
to  ry  inaccurate  drag  predictions  for  the  fully  developed  condition. 
However,  reasonable  friction  factors  can  be  obtained  if  the  modification  of 
equation  (2-46)  is  used. 

5.  Pressure  signatures  using  the  present  theory  compared  well  with  the 
experimental  results  of  reference  6.  Refer  to  figure  3-12. 

The  integral  method  presented  in  this  study  is  limited  to  a  right  cir¬ 
cular  cylinder  confined  in  a  concentric  tube.  An  investigation  using  the 
finite  element  or  finite  difference  solution  for  the  entire  body,  which 
includes  the  flow  over  the  ends  of  the  vehicle,  should  be  pursued.  Such  a 
solution  applied  to  the  confined  wake  would  be: 

1 .  More  accurate 

2.  Not  limited  to  Lv/d  >  10  (limited  by  Hoerner  coefficient) 

3.  Applicable  to  a  variety  of  tail  lengths  and  degrees  of  tapering. 


Also,  this  solution  applied  to  the  flow  over  the  nose  would  yield: 


1.  More  accurate  pressure  signatures  over  the  nose 

2.  Solutions  for  a  variety  of  nose  shapes 

3.  Much  better  approximation  of  the  velocity  profile  at  the  entrance 
of  the  annulus. 

It  is  not  certain  how  much  error  was  introduced  by  assuming  a  uniform 
velocity  profile  at  the  annular  entrance.  Intuitively,  it  is  realized  that 
the  profile  will  actually  be  skewed,  and  that  the  effect  of  skewing  will  be 
felt  throughout  the  annulus.  The  fluid's  momentum  has  not  been  properly 
accounted  for  in  this  study  since  the  boundary  layer  on  the  vehicle  surface 
was  assumed  to  begin  after  the  nose  section.  Moreover,  approximately 
one-half  of  the  pressure  drag  can  be  attributed  to  the  pressure  drop  over 
the  nose  of  the  vehicle  (figure  3-12).  The  use  of  continuity  and  Bernoulli 
equations  to  adjust  C  (nose  coefficient  from  Hoerner)  for  the  confined 
flow  problem  is  only  an  approximation  for  the  pressure  drop  across  the  nose 
(see  appendix  H).  A  more  accurate  method  for  determining  this  pressure 


field  should  be  pursued. 

The  coefficient  C,  (base  coefficient)  was  assumed  by  Hoerner  to  be 

_  *  ®  if 

a  function  of  (x^)e  and  (x^)^  as  s^own  appendix  H.  This 
relationship  may  be  true  for  external  flows  only  and,  therefore,  may  be 
inappropriate  for  internal  flows. 

The  analysis  presented  in  this  report  is  limited  to  vehicles  that  have 


reached  constant  velocities.  This  may  be  largely  the  case  for  such  appli¬ 
cations  as  trains  moving  through  tunnels.  It  is  not  exactly  the  situation 
encountered  in  a  launching  process  where  the  vehicles  or  projectiles  are 
accelerating  while  in  the  tube.  This  analysis  should  be  extended  to  include 
vehicle  accelerations. 


Entry  lengths  L  /D^,  reported  in  reference  32,  have  been  improved 
by  10-40  percent.  The  following  factors  contributed  to  this  improvement: 


WWW sew 


1.  First-order  (instead  of  second-order)  curve  fits  of  the  last  four 
computed  points  of  the  developing  region  were  used  to  determine  the  slope  of 
the  boundary  layers  "i"  and  "o"  at  the  last  computed  point,  to  extrapolate 
for  the  developed  values. 

2.  Double  precision  was  used  when  performing  the  momentum  balances, 
equations  (3-11)  through  (3-24). 
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APPENDIX  A 


COMPUTATIONAL  SCHEME  FOR  DETERMINING 
DRAG  IN  THE  REGION  OF  FULLY  DEVELOPED  FLOW 


The  velocity  profiles  were  generated  following  the  scheme  in  figure  A-l. 
The  generated  profiles  must  satisfy  continuity  and  the  imposed  boundary 
conditions  listed  in  section  2.  Velocity  profiles  were  first  generated  for 


region  "o,"  then  for  region  "i."  These  profiles  must  meet  at  r  *  r  ,  the 

location  of  maximum  velocity.  This  means  that  at  y  -  y  ,  the  velocity  U 

mo  mo 


predicted  by  region  "o"  computations  must  closely  match  the  velocity  U  at 
y+  =  y^\  ,  predicted  from  region  "i"  computations. 


The  general  iteration  procedure  was  obtained  from  reference  U.  Graphic 
representation  of  the  iterative  scheme  used,  figure  A-l,  will  now  be 
explained  in  greater  detail.  Numbers  preceding  the  following  paragraphs 
correspond  to  the  block  numbers  in  figure  A-l. 


© 


The  quantity  r*  was  assumed,  which  is  equivalent  to  specifying  the 
shear  stress  on  the  tube  wall  from  r+  *  r 


-vlx  /p  .  An  alternate  physical 
+  w  ^  M  o 

implication  of  the  choice  for  rQ  is  apparent  from  the  relationship 
t  *  -|i(du/dy),  which  implies  that  the  slope  of  the  velocity  profile  has 


been  assumed.  Having  a  good  initial  guess  for  r*  would  minimize  the  number 
of  iterations  required  for  convergence.  The  first  estimate  of  r*  for  a 


specified  Utu^e  was  found  as  described  in  the  following  paragraphs. 


The  bulk  velocity  U^n  entering  the  annulus  is  obtained  from  conti¬ 
nuity,  with  the  assumption  that  the  fluid  acts  like  slug  flow  at  the  annular 


entrance.  The  volume  flux  Q  that  approaches  the  stationary  vehicle  is 
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Figure  A-l.  Soheme  for  Determining  Vehiole  Drag 
In  Fully  Developed  Flow 
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From  continuity,  Q.  =  Q  ;  the  relationship  between  the  inlet 

1  tl  & 

velocity  U.  and  the  velocity  U.  ,  is  then 
'  in  tuoe 


Utube  b 
(b2  -  1) 


(A-l) 


For  turbulent  annular  flow  through  a  concentric  annulus,  White  (reference  32) 
suggests  the  friction  factor  is  of  the  form 


_  64£ 


(r  -r.)  (r2-r2) 
_ o  i  V.  o  l  J _ 

r4-r4-(r2-r?')  log  (r  /r.) 
o  i  V  o  i  J  Beoi 


The  shear  stress  on  the  tube  wall  is  found  from  the  definition  of  friction 


factor,  assuming  that  *  t  .  Then, 


fplTn  (b-1) 
2(b+l) 


Once  the  wall  shear  is  known,  rQ  may  be  easily  determined  from  equation 
(2-2a)  and  one  similar  to  equation  (2-2c). 

®A  value  of  a  ■  r  ./r,  is  specified  next;  r  .  is  obtained  from 

mu  1  rat 

equation  (2-14). 

©The  quantities  x  ,  x . ,  y*  ,  y+. ,  x/x  ,  x/x.  are  determined  from 

o  l  no  ni  o  l 

the  equations  specified  in  block  (3)  of  figure  A-l.  These  values  may  be 
determined  once  values  for  r*  and  a  have  been  specified.  The  quantities 
x/x ^  and  xlxQ  are  needed  before  equations  (2-8)  and  (2-9)  <*an  be  solved. 

Region  "o"  velocity  profiles  are  generated  by  integrating  equations 
(2-9a)  and  (2-9b).  The  Runge-Kutta  method  is  used  to  perform  the  integra¬ 
tion  as  explained  in  section  2.  The  step  size  for  numerical  integration  is 
determined  from  a  parametric  study  of  8  >  fcn(Rev,b).  The  results  of  the 
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study  are  shown  in  figure  2-5,  from  which  a  curve  fit  was  made.  Only  100 
steps  were  necessary  to  obtain  sufficient  accuracy  when  integrating  for  the 
sublayer  profiles.  Therefore,  the  curve  fit  applies  only  to  the  turbulent 
layer  and  will  predict  velocity  profiles  within  ♦  3.5  percent.  The  curve 

4 

fit  applies  for  the  radius  ratio  range  1.01  <  b  <  2.0  and  for  10  < 

g 

Re  <  10  ,  where  Re  is  based  on  vehicle  radius  and  speed. 


© 


5 J  Equation  ( 2— 20d )  was  derived  by  assuming  Utube  =  0.0.  This  value 
was  assumed  to  simplify  the  algebra  during  the  derivation.  At  this 
point,  the  velocity  profile  must  be  corrected  by  adding  the  dimensionless 
value  for  a  chosen  tube  velocity  t0  the  existing  velocity 

components  u*^: 


+  +  '  ..  .  * 
u„=u„+U,/u 
o2  o2  tube  o 


® 


6)  The  quantity  Umo,  the  dimensional  velocity  at  the  edge  of  the 

turbulent  boundary  layer  for  region  "o,”  must  be  larger  than  the  velocity 

CJin  entering  the  annulus.  This  is  true  because  the  boundary  layers  have 

grown  and  met,  causing  the  core  fluid  to  accelerate  to  fully  developed 

velocities  so  that  continuity  could  be  satisfied.  The  quantity  is 

checked  against  U.  to  accelerate  convergence.  If  U  <  U.  ,  the  assumed 
in  *  mo  in 

value  for  the  slope  of  the  velocity  profile  at  the  tube  wall  is  considered 
to  be  too  small;  thus  r*  is  increased.  A  good  control  on  the  amount  of 
increase  was  found  to  be 


<  ■  1A<(Uln'W  ■ 


Blocks  (3)  -  (7)  are  then  repeated  until  U  >  U. 

mo  m 

( 1 ^  Once  the  velocity  profile  for  region  "i"  is  computed 

by  integrating  equations  (2-8)  in  a  similar  manner  as  was  done  for  region 
"o"  (block  4). 
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®The  dimensional  velocities  U  and  U  .  must  be  nearly  equal  to 

mo  mi  1  n 

comply  with  boundary  condition  no.  4  in  section  2. 


rp\ 

! 

1 


10)  If  the  condition  in  this  block  is  not  satisfied,  the  value  for 
from  equation  (2-31)  is  compared  with  the  value  a,  where 


a,  »  r  . 
hg  mt 


/r.  =  jO . 5 ( r  -  r.)  +  r.l/r.  . 

l  L  oi  lj  i 


(ll)  If  a  <  a.  ,  then  r  is  assumed  to  be  too  low.  For  annular 
> — /  new  —  tig  o 

tube  flow,  r  will  lie  above  the  centerline  of  the  gap.  This  is  especially 
true  if  the  tube  is  moving.  By  increasing  r*,  tq  is  really  being 
increased,  which  implies  that  the  slope  of  the  velocity  profile  at  the  tube 
wall  is  effectively  increased.  The  quantity  r*  is  increased  by  the 

following  amount: 


& 


r+  a  r+  +  10,000  . 
o  o 


This  rate  of  increase  for  r*  proved  to  work  well  in  speeding  the  conver¬ 
gence  for  the  range  of  Reynolds  numbers  and  radius  ratios  tested.  Blocks 
(3)-(10)  are  repeated  until  block  (10)  is  satisfied. 


(12)  If  a  >  a.  ,  then  r  must  be  located  above  r  .  *  a.  r . , 
vy  new  -  hg  mt  mt  hg  i 

which  is  the  expected^  condition.  is  checked  to  see  if  it  was  within 

10  percent  of  U  .  ' 


13)  Since  is  not  within  10  percent  of  U  ,  the  full  value  of 
from  equation  (2-31)  is  taken  and  blocks  (3)-(12)  are  repeated  until 


new 

block  (12)  is  satisfied. 


(14)  Since  is  within  10  percent  of  Uno  and  block  (9)  is  not 
satisfied,  a  more  precise  value  for  a  must  be  determined.  This  is  done  by 
taking  an  average  between  a  and  ®new*  Blocks  ( 3 )— ( 9 )  are  repeated  until 
block  (9)  is  satisfied. 


i|J 


(l5)  Although  U  ^  is  within  O.i  percent  of  UmQ,  continuity  must  be 
checked.  The  voLume  flux  that  enters  the  annulus  must  be  approximately 
equal  to  the  volume  flux  at  any  axial  location  in  the  annulus.  This  is 
achieved  by  calculating  the  bulk  velocity  from  equation  (2-35). 

@  The  bulk  velocities  are  used  as  a  basis  fcr  checking  continuity 
rather  than  Reynolds  numbers  so  that  greater  sensitivity  can  be  gained  from 
the  smaller  magnitudes  of  the  velocities  compared  to  the  typically 
large  values  of  Reynolds  numbers. 

©Since  U,  is  not  within  0.1  percent  of  U.  ,  a  new  r+  (slope  of  the 
b  in  o 

velocity  profile  at  the  tube  wall)  is  calculated  from  equation  (2-37). 
Blocks  (3)-(7)  are  repeated  until  block  (17)  is  satisfied. 


(1^)  Since  the  bulk  velocity  in  the  annulus  compares  well  with  the 

bulk  velocity  entering  the  annulus  U.  ,  the  solutions  for  friction  factor 

in 

and  wall  shear  stress  have  been  found.  Therefore,  the  total  drag  may  be 
calculated  from  equation  (2-42). 


APPENDIX  B 


TABULATED  RESULTS  OF  ANALYSIS 
FOR  FULLY  DEVELOPED  FLOW 


Re 

V 

b 

a 

ffd 

fM 

(t/ 

fd 

(x  )* 

0  fd 

"fd/4P>I 

104 

1.01 

1.00505 

101. 

0.03902 

25.09 

24.51 

1.03 

103 

1.00506 

49.3 

0.02004 

12.32 

12.00 

0.97 

106 

1.00506 

29.5 

0.007713 

7.469 

7.253 

0.99 

107 

1.00506 

19.7 

0.005438 

4.942 

5.000 

1.01 

108 

1.00506 

14.0 

0.005929 

3.559 

4.848 

1.02 

104 

1.02 

1.01019 

25.5 

0.03896 

5.927 

1.03 

105 

1.01023 

12.4 

0.02002 

2.850 

0.97 

106 

1.01024 

7.43 

0.01181 

1.901 

0.99 

107 

1.01025 

4.97 

0.007707 

1.277 

1.01 

108 

1.01026 

3.54 

0.005435 

0.9110 

1.02 

104 

1.03 

11.4 

0.03889 

2.658 

1.02 

105 

5.55 

0.01999 

1.312 

0.97 

106 

3.33 

0,01180 

0.8667 

0.99 

107 

2.23 

0.007700 

0.5547 

1.01 

108 

1.59 

0.005431 

0.4146 

1.02 

104 

1.04 

1.02076 

6.51 

0.03883 

1.509 

1.03 

105 

1.02091 

3.15 

0.01997 

0.7637 

0.97 

106 

1.02095 

1.89 

0.01179 

0.4996 

0.99 

107 

1.02100 

1.25 

0.007694 

0.3353 

1.00 

108 

1.02106 

0.901 

0.005427 

0.2397 

1.02 

4.21 
2.03 

1.22 
0.817 
0.584 


0.03876 

0.01994 

0.01178 

0.007687 

0.005423 


1.090 

0.5381 

0.3275 

0.2201 

0.1574 


0.9730 

0.4713 

0.2838 

0.1892 

0.1347 
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p» 
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104 

1.06 

1.03171 

2.95 

0.03870 

0.6867 

1.03 

•s 

103 

1.03199 

1.42 

0.01991 

0.3630 

0.97 

r 

LO6 

1.03212 

0.856 

0.01176 

0.2304 

0.99 

LO7 

1.03226 

0.573 

0.007681 

0.1566 

1.01 

108 

1.03238 

0.408 

0.005419 

0.1104 

1.02 

•  ? 

104 

1.07 

2.19 

0.03863 

0.5126 

1.03 

► 

i 

LO3 

1.05 

0.01989 

0.2756 

0.96 

t 

106 

0.634 

0.01175 

0.1712 

0.98 

107 

0.429 

0.007674 

0.1175 

1.02 

104 

1.08 

1.04295 

1.69 

0.03857 

0.4051 

1.03 

l 

105 

1.04345 

0.811 

0.01986 

0.2179 

0.96 

106 

1.04384 

0.490 

0.01174 

0.1326 

0.98 

107 

1.04398 

0.328 

0.007667 

0.090215 

1.01 

LO8 

1.04445 

0.234 

0.005412 

0.06010 

1.02 

« 

\ 

\ 

f 

LO4 

1.09 

1.35 

0.03850 

0.3228 

1.04 

LO5 

0.645 

0.01984 

0.1378 

0.96 

*  h 

106 

0.391 

0.01173 

0.1061 

0.99 

107 

0.262 

0.007661 

0.07453 

1.01 

108 

0.187 

0.005408 

0.05338 

1.02 

. 

104 

1.1 

1.1 

0.03844 

0.02757 

1.03 

105 

0.529 

0.01981 

0.1493 

0.96 

LO6 

0.320 

0.01171 

0.08707 

0.99 

107 

0.215 

0.007654 

0.06096 

1.01 

108 

0.153 

0.005404 

0.04433 

1.02 

i 

104 

1.07999 

0.05497 

0.03801 

0.1562 

0.1151 

1.02 

105 

1.143 

1.08199 

0.2677 

0.01964 

0.07994 

0.05561 

0.96 

106 

1.08303 

0.1623 

0.01163 

0.04902 

0.03310 

0.99 

107 

1.08381 

0.1089 

0.007610 

0.03319 

0.02191 

1.01 

108 

1.08433 

0.02399 

0.0156 

1.36E5 

1.164 

1.0974 

0.1872 

0.01822 

0.05722 

0.03747 

0.93 
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LO4 

105 

106 

107 

LO8 


0.2994 

0.03779 

0.0824 

1.03 

0.144 

0.01955 

0.04584 

0.96 

0.0870 

0.01159 

0.02594 

0.98 

0.0588 

0.007588 

0.01911 

1.01 

0.0419 

0.005366 

0.01380 

1.02 

104 

1.333 

1.2075 

0.1206 

0.03694 

0.04121 

0.02108 

1.03 

105 

1.21725 

0.0583 

0.01920 

0.02108 

0.009693 

0.95 

106 

1.2221 

0.0351 

0.01142 

0.01311 

0.00565 

0.96 

107 

1.22576 

0.0237 

0.007500 

0.00903 

0.003699 

0.99 

108 

0.0172 

0.005314 

0.006542 

1.01 

104 

1.4 

0.0886 

0.03652 

0.03170 

1.02 

105 

0.01628 

0.01903 

0.01642 

0.95 

106 

0.0260 

0.01134 

0.01021 

0.96 

107 

0.0177 

0.007456 

0.007128 

1.00 

108 

0.0126 

0.005289 

0.004927 

1.00 

104 

1.6 

1.41288 

0.0464 

0.03532 

0.01929 

0.006823 

1.01 

105 

1.43954 

0.0223 

0.01854 

0.009974 

0.002835 

0.92 

106 

1.45194 

0.0137 

0. 01111 

0.006299 

0.001606 

0.95 

107 

1.46149 

0.00921 

0.007329 

0.004430 

0.001034 

0.97 

108 

0.00667 

0.005215 

0.003038 

0.98 

104 

1.8 

0.0304 

0.03422 

0.01268 

0.99 

105 

0.0146 

0.01809 

0.007386 

0.90 

106 

0.00896 

0.01089 

0.004720 

0.92 

107 

0.00620 

0.007211 

0.003315 

0.96 

LO8 

. 

0.00440 

0.005146 

0.002370 

0.95 

104 

2.0 

0.0223 

0.03322 

0.01122 

0.00255 

0.97 

105 

1.80723 

0.0107 

0.01767 

0.006010 

0.0009731 

0.87 

106 

i. 83229 

0.00662 

0.01068 

0.003892 

0.0005305 

0.89 

107 

1.84989 

0.00452 

0.007100 

0.002794 

0.0003334 

0.92 

108 

0.00332 

0.005081 

0.001988 

0.94 
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APPENDIX  C 


COMPARISON  OF  MOODY  CHART  METHOD  AND 
PRESENT  THEORY  FOR  PREDICTING  PRESSURE  DROP 
THROUGH  ANNULI  WITH  A  MOVING  BOUNDARY 


The  following  example  was  performed  to  illustrate  the  difference  between 
the  pressure  drop  from  the  fully  developed  analysis  and  that  from  using  the 
Moody  chart. 

The  comparison  applies  for  the  following  condition: 

Vehicle  roughness  =  smooth, 

r.  =  0.875  ft, 

1 

b  =  1.1429, 

L  -  21  ft, 
v 

Uv  ■  40  ft/s, 

Fluid  *  sea  water  @  60°F,  v  *  0.1258  x  10  ^  ft2/s. 

Then, 

Rev  =  Uvr./v  *  2.782  x  106  .  (C-l) 

The  Reynolds  number  associated  with  the  Moody  chart  is 

Reb  *  UfeDh/v  -  3.39  x  106  .  (C-2) 

By  definition,  the  actual  fully  developed  pressure  drop  is  related  to  the 
friction  factor  f^  by 


Ap 


fd 


P^vL 


fd 


2*cDh 


(C-3) 


In  comparison,  a  computation  using  the  Moody  chart  would  result  in  a 
pressure  drop  given  by 


where,  neglecting  entrance  and  exit  losses,  the  Moody  head  loss  is  related 
to  an  appropriate  average  kinetic  energy  in  the  passage: 


2 

,  fMLUav 

£  '  Dk2g 


(C-4) 


From  the  continuity  relation,  the  bulk  and  vehicle  velocities  are  related  by 


Uu  »  U  b2/(b2  -  1)  . 

D  V 


Because  the  inner  and  outer  tube  walls  are  moving  at  different  speeds,  the 
inner  shear  stress  correlates  with  U.  ,  and  the  outer  shear  stress 

D 

correlates  with  (U^  -  Uv).  By  choosing  the  arithmetic  mean  of  these  as 
the  appropriate  average  velocity  to  use  in  equation  (C-4), 


U 

av 


V  2<b2-l> 


(C-5 ) 


Combining  equations  (C-4)  and  (C-5)  and  using  the  Moody  chart  yields 
the  pressure  drop: 


AP, 


M1 


i  P^L 


M 


2*cDh 


2 


\  +1 
2(b2-l) 


(C-6 ) 


where  is  computed  from  equation  (2-46)  using  Re^  as  the  appropriate 
Reynolds  number.  The  two  estimates  of  pressure  drop  given  by  equations  (C— 3 ) 
and  (C-6)  should  be  equal  if  the  assumptions  in  the  comparison  are 
realistic.  The  ratio  of  these  two  equations  is 


For  Che  particular  example  given  in  this  appendix,  the  author  has  computed 

,6 


=  0.L35.  From  equation  (2-46),  with  Re^  =  3.39  x 
With  b  =  1.143,  equation  (C-7)  predicts 


10 


fM  =  0.00933. 


Ap 


fd 


^Pv 


0.135 

0.00933 


(0.2657)4  =  1.02 


For  this  example,  the  error  in  the  Moody  estimate  is  only  2  percent.  This 
ratio,  tabulated  in  table  B-l  for  a  variety  of  values  of  Re  and  b,  shows 
that  the  Moody  approximation  will  result  in  an  average  error  of  +  1.3 
percent.  The  maximum  error  of  13  percent  occurred  at  b  *  2.0,  Re^  =  10^. 
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APPENDIX  D 


COMPUTATIONAL  SCHEME  FOR  DETERMINING 
DRAG  IN  THE  REGION  OF  DEVELOPING  FLOW 


The  scheme  used  in  this  study  for  determining  drag  in  the  region  of 
developing  flow  is  depicted  in  figure  D-l.  The  following  text  provides 
detailed  explanations  of  the  scheme.  (Numbers  preceding  paragraphs 
correspond  to  the  block  numbers  in  figure  D-l.) 

(T)  These  inputs  are  the  only  ones  necessary  to  run  the  computer  code. 
NUMX  is  the  number  of  axial  stations  to  be  used  in  the  analysis.  A  typical 
value  was  NUMX  =  20.  However,  greater  accuracy  may  be  obtained  by  assigning 
larger  values  to  NUMX  at  the  expense  of  increased  central  processing  unit 
(CPU)  time.  The  longest  CPU  time  experienced  by  the  author  for  a  solution 
to  the  developing  velocity  profiles  was  approximately  2  hours.  This  was 
for  just  20  axial  stations,  and  Rev  *  10^  for  b  »  2.0. 


© 


3  2 

The  fluid  density  p  (lbm/ft  )  and  kinematic  viscosity  v  (ft  /s) 


were  acquired  from  subroutine  "Viscosity."  Sutherland's  law  was  used  to 
compute  p  for  air,  and  the  ideal  gas  equation  was  used  to  find  p. 

The  viscosity  for  water  was  found  by  using  an  empirical  formula  developed  by 
Bingham  from  reference  18. 
from  equations  (2-43)  and  (2-44). 


The  values  for  a  and  S.  were  obtained 

lmx 

If  Rev  is  something  other  than  the  values 


for  which  the  equations  were  fitted,  a  linear  fit  between  the  encompassing 
Reynolds  numbers  will-  provide  a  fair  estimate  for  6^^.  The  number  of 
iterations  used  in  the  Runge-Kutta  routine  for  integrating  the  sublayer 
equations  (2-8a)  and  (2-9a)  was  set  to  ITER  ■  1000.  This  was  found  to  be 
sufficiently  large  for  all  radius  ratios  and  Reynolds  numbers  tested.  The 
number  of  iterations  used  to  integrate  the  turbulent  layer  equations,  ITER2, 
was  determined  from  the  curve  fit  shown  in  figure  2-5.  The  limit  of  integra¬ 
tion  for  the  region  "i"  turbulent  layer,  for  each  station,  depends 
on  the  number  of  stations  NUMX  to  be  evaluated.  For  example,  if  10  stations 
were  to  be  used,  6^  would  be  divided  into  10  boundary  layer  heights 
(figure  D-2).  The  first  height  would  be  fit,  *  &,  /NUMX.  The  next 

+  -4s 

height  would  be  ■  Sjj,  +  8 ^  /NUMX,  and  so  on.  Each  height  corresponds 


D-l 


CALL  (UMOUTtM  "IMT1Q  T 
TO  INTMMATI  CO.  (M) 
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Figure  D-l.  Soheme  for  Determining  Vehiole  Drag 
in  Developing  Flow 


Figure  D-l.  Seheme  for  Determining  Vehlole  Drag 
In  Developing  Flow  (Coni’ d) 


*1 

I 


to  a  unique  axial  location.  For  each  S+  ,  values  of  S+  and  r+  were 
n  moo 

estimated  using  the  following  curve  fits  of  data  from  reference  5: 


S  +  =  0.41316T  , 
o  1 


r+  =  ABS(8/9rt  -  1333.33)  . 
o  1 


(D-l ) 


(D-2) 


The  values  obtained  from  these  curve  fits  serve  only  as  first  estimates  and 
will  be  refined  by  iteration. 

Subroutine  "INTEG"  at  this  point  integrates  equation  (2-8a)  from 
0  to  yt^.  The  Runge-Kutta  method  is  used  to  numerically  integrate  the 
equation.  Subroutine  "INTEG"  outputs  the  following:  u^  (y^.),  ut^  (y*^)f 
and  dut^/dyT^.  The  quantities  ut^,  yt^  and  dut^/dyt^  are  used  as 
inputs  to  subroutine  "INTEG2,"  which  integrates  equation  ( 2— 8b ) . 

(T)  Subroutine  "INTEG2"  is  called  at  this  point  to  numerically  integrate 
equation  (2-8b)  from  yt^  to  stn«  The  outputs  of  this  routine  are: 

tt2i<!r2i>*  “mi  *  u2i^imt 

(T)  The  dimensional  velocity  U,^  at  the  boundary  layer  edge  of  region 
"i"  is  determined  from 


+  * 

U..  *  u  .  u. 

Si  mi  i  - 


(T)  The  following  condition  is  used  to  accelerate  convergence: 


Ut.  >  U.  . 
Si  in 


The  flow  must  accelerate  through  the  annular  core  since  the  fluid  is  slowed 
down  within  the  growing  viscous  boundary  layers  emanating  from  both  the 
vehicle  and  tube  walls. 
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(T)  Since  the  condition  in  block  (6)  has  not  been  satisfied,  it  is 


assumed  that  rt  is  too  low.  Therefore,  rt  must  be  increased  until  block  (6) 


is  satisfied,  which  effectively  means  the  wall  shear  stress  is  actually 


being  increased. 


A  new  x .  is  computed  from  equations  similar  to  equation  (2-2a). 
Blocks  (5)  through  (7)  are  repeated  until  block  (6)  is  satisfied. 


(^8^  Since  block  (6)  has  been  satisfied,  the  velocity  at  the  present 
location  must  be  greater  than  the  values  at  the  previous  locations.  As  the 
code  progresses  into  the  annulus,  as  shown  in  figure  (D-2),  values  of 


6.  will  increase, 
in 


(^9 )  Since  the  condition  in  block  (8)  has  not  been  satisfied,  rt  is 
increased  by  0.001  percent.  Blocks  (5)  through  (8)  must  be  repeated  until 


this  condition  is  satisfied. 


UJ)  Now  that  a  reasonably  large  value  of  Ufii  has  been  found,  the 


velocity  profile  in  the  sublayer  for  region  "o"  is  found  by  using  subroutine 
"INTEG."  This  routine  integrates  equation  (2-9a)  between  the  limits  of  0  to 


y*,-  "INTEG"  yields:  u*o(y*Q),  u*(y+,),  and  (du*/dy*)„. 


(Ty  The  values  u^(y^)  and  (du*/dy*)^  are  used  as  inputs  to 
subroutine  "INTEG2,"  which  numerically  integrates  equation  (2-9b)  between 
the  limit,  of  y0%  to  S^.  "INTEG2"  yield,  [«^.  y*0>  -  <(y^t)]  - 


12)  The  assumed  tube  velocity  U  is  nondimensionalized  as  follows: 


.  /u  . 
tube  o 


( D— 3 ) 


1J)  UgQ  is  adjusted  by  to  account  for  the  moving  tube  as  follows: 


„+  +  „+ 

U.  ■  u  +  U.  . 
So  mo  tube 


iffl 

ty 

$ 


w 


It  is  not  necessary  to  adjust  the  whole  region  "o"  profile  at  this  point 

since  more  iterations  may  be  required  to  match  U.  and  U... 

So  Si 


and  must  be  nearly  equal  due  to  the  inviscid  nature  of  the 
potential  core.  The  fluid  in  the  core  region  has  not  experienced  any  viscous 
effects  due  to  the  presence  of  the  walls.  The  term 


CHK1  =■ 


USi  -  USo 
USi 


provides  a  check  on  how  close  Uf  is  to  U...  It  also  indicates  the  direction 

60  61 

of  correction.  The  matching  criteria  of  0.001  percent  is  chosen  to  ensure 
numerical  stability  during  the  computations  for  Refc,  equation  (2-34),  and 
dp/dx,  equation  (3-25). 


(l6)  A  series  of  checks  is  performed  on  the  value  of  r+.  Based  on  the 

.  +  .  0 
disagreement  between  and  UgQ,  rQ  is  adjusted  by  using  a  linear  interpo¬ 
lation  scheme.  Blocks  (12-16)  are  repeated  until  block  (15)  is  satisfied. 

©  At  this  point,  block  (15)  has  been  satisfied,  and  now  all  the  veloc 
ities  in  the  boundary  layer  of  region  "o"  must  be  scaled  by  the  dimension¬ 
less  tube  velocity  u£ube  to  account  for  the  moving  tube. 

^8)  At  this  point,  the  core  velocities  and  UgQ  are  nearly  equal 
(within  0.001  percent).  Continuity  must  now  be  checked  by  computing 
equation  (2-34).  Barts  1  and  2  of  that  equation  are  integrated  using  the 
trapezoidal  rule.  The  bulk  velocity  is  used  rather  than  the  Reynolds  number 
to  check  for  continuity.  Bulk  velocities  are  generally  much  smaller  in 
magnitude  than  the  typically  large  values  of  Reynolds  numbers,  which  may 
make  the  check  too  sensitive  to  differences  of  large  numbers. 

©)  If  continuity  is  satisfied,  the  location  of  velocity  profiles  is 
predicted.  If  continuity  is  not  satisfied,  then  r£  must  be  adjusted  until 
continuity  is  satisfied. 


(20)  The  quantity  is  adjusted  by  using  a  linear  interpolation  scheme 
based  on  the  resulting  value  of  from  the  previous  guess  of  rt.  Increas¬ 
ing  rt  will  increase  U,  . 

1  b 


(2^  The  first  velocity  profile  within  region  "i"  is  located  by  solving 
equations  (3-11)  and  (3-16)  simultaneously.  Successive  profiles  within  this 
region  are  located  by  solving  equations  (3-23)  and  (3-16)  simultaneously. 


(22)  The  first  velocity  profile  in  region  "0"  is  located  by  solving 
equations  (3-12)  and  (3-16)  simultaneously.  Successive  profiles  within  this 
region  are  located  by  solving  equations  (3-24)  and  (3-16)  simultaneously. 


(23)  Now  the  velocity  profile  for  region  "i’’  is  checked  to  see  if  it  has 
the  same  location  as  that  for  region  "0." 


24)  The  quantity  Sq  is  adjusted  based  on  the  discrepancy  between 


(x)^  and  (x)q  and  the  previous  guess  of  $*. 


(25)  Now  both  t^(x)  and  Ap(x)  can  be  integrated  over  the  entry  length  L  . 
These  values  are  then  multiplied  by  the  vehicle  surface  area  within  the  devel¬ 
oping  region  and  the  projected  vehicle  area,  respectively,  to  obtain  the  drag 
shown  by  equation  (3-26). 
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APPENDIX  E 

DRAG  PREDICTION  USING  FULLY  DEVELOPED  AND  DEVELOPING 
FLOW  RESULTS  TO  QUANTIFY  THE  VARIOUS  DRAG  CONTRIBUTIONS 


This  appendix  demonstrates  how  one  would  compute  the  drag  from  both  the 
developing  and  fully  developed  flows  through  the  annular  gap.  This  is  accom¬ 
plished  by  using  only  the  graphs  generated  as  a  result  of  this  investigation. 
Consider  the  following  case: 


r.  =  0.875  ft, 

1 

b  =  1.143, 

L  *  21  ft, 

V 

Uv  =  40  ft/s. 

Fluid  =  sea  water  @  60°F,  v  a  0.1258  x  10  ^  ft2/s. 
Re  a  U  r./v  =  2.782  x  106  . 

V  VI 


The  total  drag  force  on  the  vehicle  is: 
0 


(DRG)t  +  (PDRG)@  +  (SHRDRG)e  +  (PDRG)fd 


+  (SHRDRG) 


fd 


‘X 


(E-l ) 


Consider  only  the  drag  contribution  due  to  the  flow  through  the  annulus. 
The  nose  drag  (F^)  and  base  drag  (Fg)  are  not  included  in  this 
discussion.  They  can  be  approximated  as  described  in  appendix  H. 

The  drag  forces  contributed  by  the  developing  portion  of  the  flow  are: 


*  pU2A 

(PDRG)  -  f  -r~  »  14,788  lbf  , 

e  e  Dh  2gc 


E-l 


where  f  =  0.1646  (figure  3-19),  and  L  /D^  =  23.46  (figure  3-22),  and 


_  *  puv 

(SHRDRG)  =  (t.)  r —  2irr .  L  =  2354  Ibf  , 
e  l  e  2g=  l 


where  (x£)g  *  0.04584  (figure  3-20). 


The  drag  forces  contributed  by  the  fully  developed  portion  of  the  flow 


(L  -L*)  p(J2A 

(tDRG)fd  *  -V-  Efd  -if  =  31-252  lb£  • 

h  c 


where  f£d  =  0.1348,  and 


*  PV  * 

( SHRDRG )fd  =  (Ti)£d  2irr.(Lv  -  L  )  »  5424  lbf  , 


where  C*c^ ) =  0.04094  (figure  2-10).  Then,  from  equation  (E-l), 


(DRG)t  *  53,818  lbf. 


For  the  case  presented  here,  the  drag  contributions  can  be  summarized 


as  follows: 


(PDRG)e  »  22.5  percent 


( SHRDRG )e  *  4.4  percent 


Drag  due  to  developing  flow. 


(PDRG)fd  ■  58.1  percent 


(SHRDRG)fd  *  10.1  percent 


Drag  due  to  fully  developed  flow. 


The  shear  and  form  drag  contributions  for  developing  and  fully  developed 


flows  will  vary  with  Lv,  b,  and  Rev.  The  shear  drag  will  usually  be 


less  significant  than  the  form  drag  for  both  flow  regions. 
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APPENDIX  F 

TABULATED  RESULTS  OF  ANALYSIS 
FOR  DEVELOPING  FLOW 
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APPENDIX  G 


DATA  TABLES  FOR  CJV  =  40  FEET/SECOND,  b  =  1. 143  RL'N 


Table  G 

-1.  Entry  Region  Result 

.  for  U 

V 

=  40  ft/s, 

r.  =  10.5 
i 

f! 

O 

U 

c 

•r( 

12.0  in 
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s. 

1 

$ 
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Ti 

T 
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dp/dx 

f 

e 

(ft) 

(ft/s) 

(ft) 

(ft) 

(lbf/ft2) 

(lbf/ft2) 

(lbf/ft3) 

0.000 

0.0000 

170.5273 

0.000000 

0.000000 

370.339569 

248.992249 

0.067 

0.2694 

171.0750 

0.002395 
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109.086792 
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1.751 

6.9968 

176.2685 

0.032868 
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0.036192 
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0.053309 
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68.722549 
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0.043914 

68.371521 
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-928.36432 

0.1459 
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23.4555 

185.0870 

0.071477 

0.053648 

66.070000 

43.599998 

-858.90002 

0.1350 

Integrated  Average 

72.99 

47.57 

-1048. 

0 . 1646 

Fully 

Developed — >  182.75  0.07287  0.05213 

Analysis 

65.197 

43.599 

-858.87 

0.1348 

XDifference — >  1.3  1.9  2.9 
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Developing  and 

Fully  Developed 

1.3 

0.0 
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Table  G; 

-2.  Volume  Flux  for  Profiles 

Shown 

in  Figure  3-11 

Profile 

# 

*1 

^CORE 

*0 

QTOt 

(ft3/s) 

<ft3/s) 

Z  of 
Total 

(ft3/s) 

Z  of 
Total 

(ft3/s) 

Z  of 
Total 

(ft3/s) 

Z  of 
Total 

1 

125.66 

5.93 

4.7 

113.40 

90.2 

5.32 

4.2 

124.65 

99.2 

2 

125.66 

36.05 

28.7 

57.41 

45.7 

32.15 

25.6 

125.62 

100.0 

3 

125.66 

61.26 

48.3 

9.34 

7.4 

54.02 

43.0 

125.62 

99.2 
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APPENDIX  t: 


METHOD  USED  TO  ESTIMATE  NOSE  AND  BASE  DRAG 


The  Bernoulli  equation  for  constant  elevation  is 


(u2  "  ul}  =  "(p2  ~  pi> 


( H—  1 ) 


where  station  1  designates  a  point  in  the  flow  just  ahead  of  the  vehicle. 
Station  2  designates  a  point  within  the  annulus  at  a  distance  from  the  throat 
equal  to  the  length  of  the  nose  section. 

Continuity  requires  that 


^  ^  =  b2/(b2  -  1)  , 

U1  A2 


2  2  2 
where  A,  ■  irr  and  A„  *  ir(r  -  r.) 

1  o  2  oi 


(H-2) 


The  analysis  assumes  that  u^  =  Uv.  Substituting  equation  (H-2)  into 
(H-l)  leads  to  the  expression  for  the  pressure  drop  across  the  nose: 


4"N-Cnt^‘[(b2/<b2-l))2-1]  ' 


The  quantity  C  is  the  nose  coefficient  from  Hoerner  (reference  30,  page 
3-12)  that  accounts  for  various  nose  types.  In  general  C  can  vary 
between  -0.05  (ogive  nose)  and  1.0  (hollow  nose).  The  drag  force  due  to  the 
pressure  drop  necessary  to  accelerate  the  fluid  over  the  nose  and  into  the 
annular  throat  is 


FM  -  C  . 
N  nt 


Ay  f(b2/(b2  -  l))2  -1  . 
*e  L  -J 
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The  drag  force  due  to  the  low  pressure  region  immediately  behind  the 
vehicle  (from  Hoerner)  can  be  approximated  by 


F3=Cb2i;Av  ’ 


where  C,  =0.29 
□ 


_  *  * 

2  (t. )  L  +(t 
Lie 


iwi 


Keep  in  mind  that  the  relationships  just  derived  for  the  drag  over  the 
nose  and  base  of  the  vehicle  are  only  approximate.  Quantities  C  and 
were  obtained  from  experiments  performed  in  an  infinite  flow  field. 
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